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INTRODUCTION 


Radiation  pattern  analysis  of  on-aircraft  antennas  at  high  fre¬ 
quencies  is  the  object  of  this  research.  It  is  a  basic  study  of 
aircraft-antenna  pattern  problems  in  which  the  antenna  is  mounted 
on  the  fuselage  near  the  too  or  bottom.  Since  it  is  a  study  of 
general-type  aircraft,  the  analytic  aircraft  is  modelled  in  its  most 
basic  form.  The  fuselage  is  assumed  to  be  a  perfectly  conducting 
convex  surface.  Thus,  the  effects  the  cock-pit  and  radome  are 
neglected  at  present.  However,  based  on  the  advances  made  on  metallic 
radomes,  the  raaomes  used  in  future  aircrafts  may  to  a  good  approxi¬ 
mation  be  perfectly  conducting.  The  wings  and  horizontal  stabilizers 
are  modelled  by  ''n"  sided  flat  plates  which  lie  in  a  plane  that  is 
parallel  to  the  fuselage  axis.  The  engines  are  approximated  by  finite 
circular  cylinders.  This  is  a  rather  crude  approximation;  yet,  it  is 
shown  that  the  engines  have  little  effect  on  the  overall  pattern. 

Thus,  it  is  not  necessary  at  present  to  consider  a  more  general 
engine  model. 

The  need  for  this  type  of  solution  is  basically  two  fold.  First, 
there  are  upwards  to  200  antennas  mounted  on  a  single  aircraft.  If 
these  antennas  can  be  located  on  the  aircraft  at  the  design  stage, 
then  one  can  expect  better  performance  in  that  optimum  locations  and 
necessary  structural  changes  can  be  anticipated.  Secondly,  antenna 
systems  are  normally  added  or  changed  in  the  course  of  an  aircraft's 
useful  lifetime.  Such  relocation  or  addition  of  antennas  has  always 
required  a  great  deal  of  engineering  time  and  money.  For  example, 
it  is  not  uncommon  for  one  to  spend  six  months  building  a  model  and 
a  second  six  months  measuring  radiation  patterns  for  antennas  mounted 
at  various  locations  around  the  structure.  On  the  other  hand,  it  is 
not  inconceivable  that  one  could  accomplish  the  same  result  in  a 
fraction  of  the  time  {perhaps  a  day)  using  computer  simulated  models 
of  the  aircraft.  Once  an  optimum  region  is  determined,  the  antenna 
can  simply  oe  flight  tested  to  ascertain  its  actual  performance.  Not 
only  can  these  computer  simulated  results  be  used  to  determine  the 
location,  but  they  can,  also,  determine  the  optimum  antenna  design  for 
a  given  application.  Note  that  these  analyses  consider  the  three 
mutually  orthogonal  delta  function  sources  which  can  be  used  to  solve 
for  the  pattern  of  an  arbitrary  fuselage  mounted  antenna  simply  by 
integrating  over  the  equivalent  aperture  currents. 
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One  of  the  first  solutions  used  to  compute  on-aircraft  antenna 
patterns  were  the  modal  solutions  for  infinitely  long  circular"l,2. 
and  elliptical  13"  cylinders.  The  fuselage  was  modelled  by  a  cylinder 
whose  elliptical  cross-section  approximated  the  fuselage  cross-section 
at  the  antenna  location.  Arbitrary  antennas  were  considered  and  the 
antenna  could  be  mounted  on  or  above  the  fuselage.  Results  were  quite 
adequate  provided  the  aircraft  structure  was  not  illuminated  too 
strongly.  In  fact,  these  solutions  have  been  the  primary  high  fre¬ 
quency  analysis  to  date.  However,  with  the  desire  to  improve  system 
performance,  versatility,  and  coverage  the  antenn?  pattern  must  be 
shaped  for  the  desired  application  in  such  a  way  that  it  can  actually 
illuminate  the  structure  quite  strongly.  In  fact,  in  many  cases  the 
system's  performance  is  dependent  on  the  pattern  effects  of  the 
secondary  contributors.  For  example,  too  strong  a  ripple  in  a  pattern 
may  cause  a  system  to  function  on  a  secondary  lobe  rather  than  the 
desired  main  beam.  Consequently,  the  demands  of  modern  sophisticated 
applications  require  that  the  analytic  model  represent  the  actual 
aircraft  in  more  detail  than  assumed  by  the  modal  solutions  alone. 

With  the  advent  of  modern  digital  compu:ers,  one  has  been  able  to 
obtain  integral  equation  solutions  for  antenna  patterns  and  impedance 
using  moment  methods.  Using  this  approach  the  surface  currents  are 
assumed  to  be  of  a  given  form  which  in  turn  radiate  in  all  directions. 
The  values  of  the  complex  currents  are  found  by  forcing  the  tangential 
component  of  the  electric  field  to  vanish  on  the  surface.  One  of  the 
first  moment  solutions  applied  to  aircraft  antenna  problems  used  the 
wire  grid  technique  which  incorporates  a  point  matching  scheme[4]. 

This  solution  requires  approximately  100  unknown  currents  per  square 
wavelength  to  be  found  in  order  that  the  wire  grid  adequately  model 
a  perfectly  conducting  surface.  A  more  sophisticated  approach  has 
been  developed  by  Richmond[5]  which  uses  a  reaction  technique  to  solve 
for  the  unknown  currents.  Yet  this  solution  still  requires  the  wire 
grid  model  of  the  aircraft  with  approximately  100  unknown  currents  per 
square  wavelength.  An  exhaustive  study  of  this  approach  has  been 
made  by  Lin[6]  in  which  he  actually  treated  the  bi static  scattering 
problem;  however,  the  same  conclusions  apply  for  the  antenna  problem., 

A  third  approach  might  be  to  divide  the  surface  into  patches  with  each 
patch  having  two  unknown  complex  currents.  Using  this  approach  per¬ 
haps  20  unknown  currents  per  square  wavelength  need  to  be  found. 
However,  all  of  these  solutions  are  restricted  to  lower  frequencies 
based  on  the  fundamental  limitation  on  the  size  of  matrices  which 
modern  computers  can  invert  without  excessive  loss  of  accuracy. 

Another  approach  that  has  found  great  succ  ss  at  solving  this 
type  of  problem  is  the  Geometrical  Theory  of  Diffraction  (GTD).  GTD 
is  basically  a  high  frequency  solution  which  is  divided  into  two 
basic  problems;  these  being  wedge  diffraction  and  curved  surface  dif¬ 
fraction.  The  wedge  diffraction  solution  has  been  applied  to  de¬ 
termine  the  radiation  patterns  of  such  basic  antennas  as  parallel 
plate  antennas[7,8,9],  parallel  plate  arrays[10,ll]  horns 
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antennas[12,13l,  parabolic  reflectorsri4,l  51,  and  rectangular  wave¬ 
guide  antennas[16l.  Both  these  diffraction  solutions  have  been 
applied  in  computing  the  pa:terns  of  antennas  mounted  on  cylin- 
ders[17, 18,19.,  rocketsr.20",  and  wingsr21".  The  only  limitation  of 
these  solutions  is  that  the  source  and  various  scattering  centers  be 
separated  by  at  least  a  wavelength.  In  some  cases  even  this  require¬ 
ment  can  be  relaxed.  Using  this  approach  one  applies  a  ray  optics 
technique  to  determine  the  fields  incident  on  the  various  scetterers. 
The  fields  diffracted  are  found  using  the  GTD  solutions  in  terms  of 
rays  which  are  summed  with  the  geometrical  optics  terms  in  the  far 
field.  The  rays  from  a  given  scatterer  tenc  to  interact  with  the 
other  structures  causing  various  high-order  terns.  In  this  way  one 
can  trace  out  the  various  possible  combinations  of  rays  that  interact 
between  scatterers  and  determine  and  include  only  the  dominant 
terms.  Thus,  one  need  only  be  concerned  with  the  important  structural 
scattering  components  and  neglect  all  other  higher-order  terms.  This 
makes  the  GTD  approach  ideal  for  a  general  high  frequency  study  of 
on-aircraft  antennas  in  that  only  the  most  basic  structural  features 
of  the  aircraft  need  to  be  modelled. 

The  basic  approach  applied  here  is  to  break  the  aircraft  up  into 
its  simplest  structural  forms.  Analyze  these  structures  using  ray 
optics  techniques  with  numerical  values  obtained  using  modal  solutions, 
physical  optics,  and  GTD.  Once  the  scattering  from  these  structures 
is  found  and  verified  by  measured  data,  they  are  adapted  to  the 
aircraft  model  simply  by  adjusting  the  incident  field.  Tn  this  way  the 
aircraft  begins  to  grow  out  of  simple  forms  into  a  structure  that 
actually  resembles  a  modern  aircraft  in  a  general  way. 

Using  these  theoretical  models  the  radiation  patterns  of  arbitrary 
antennas  mounted  on  the  fuselage  of  various  aircraft  shapes  are  com¬ 
puted.  The  principal  plane  patterns  of  Chapter  IV  include  only  those 
structural  components  whose  scattering  was  predetermined  to  be  of 
some  consequence  in  the  resulting  pattern.  In  the  volumetric  pattern 
study  of  Chapter  V  only  the  fuselage,  which  is  approximated  by  an 
arbitrary  convex  surface  of  revolution,  is  considered.  However,  based 
on  the  principal  plane  studies  the  fuselage  shape  plays  the  dominant 
effect  on  the  resulting  radiation  pattern  with  the  wings  and/or 
horizontal  stabilizers  being  the  strongest  secondary  contributors.  As 
a  result  of  these  studies,  various  near  field  scattering  problems  and 
associated  numerical  techniques  have  been  developed  which  may  find 
application  elsewhere. 
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THEORETICAL  BACKGROUND 


A.  Introduction 


As  stated  in  the  introduction  the  modal  solutions  for  the 
radiation  patterns  of  an  antenna  mounted  on  or  near  a  fuselage  hava 
previously  been  the  primary  tool  for  high  frequency  analyses  of  on- 
aircraft  antenna  problems.  In  fact,  these  solutions  can  be  quite 
adequate  provided  the  source  does  not  illuminate  much  of  the  aircraft 
structure.  Consequently,  the  approach  of  this  research  is  not  to 
replace  the  modal  solutions  for  fuselage  mounted  antenna  problems 
but  to  improve  the  solutions  by  adding  modelled  aircraft  scattering 
structures  in  the  analysis.  This  is  accomplished  by  using  the  ray 
optics  technique  which  allows  one  to  apply  various  high  frequency 
solutions  in  its  format  by  casting  these  solutions  into  ray  form. 

The  modal  solutions  for  a  plane  wave  field  incident  on  a  per¬ 
fectly  conducting  circular  cylinder  are  presented  and  used,  fre¬ 
quently,  for  numerical  computations.  The  ray  optics  techniques  will 
provide  insight  into  the  actual  mechanisms  involved  as  the  rays 
interact  with  the  cylindrical  structures.  Wedge  and  curved  surface 
diffraction  solutions,  which  constitute  the  Geometrical  Theory  of 
Diffraction  (GTD),  are  presented  in  that  they  provide  the  solutions 
necessary  to  introduce  the  additional  structures  needed  to  improve 
the  analyses.  Each  of  these  solutions  are  presented  in  basic  terms 
in  this  chapter  and  applied  to  specific  structural  scattering  problems 
in  the  following  chapters. 

B.  Modal  Solutions  for  Infinite 
Circular  Cylinders 

The  modal  solutions  for  the  total  field  at  any  point  in  space 
for  a  plane  wave  field  incident  on  a  perfectly  conducting  infinitely 
long  circular  cylinder  are  presented  in  Appendix  I.  These  solutions 
are  divided  into  the  problems  of  perpendicular  and  parallel  polari¬ 
zations  incident  on  the  cylinder  as  illustrated  in  Fig.  1.  For 
the  perpendicular  polarization  case  the  total  field  at  any  point  is 
given  by 
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where 


1  n  =  0 

2  otherwise. 


For  the  parallel  polarization  case  the  total  field  is  given  by 
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+  FIELD 


Fig.  1— Plane  wave  incident  on  an  infinitely  long 
cylinder  of  radius  (a). 


The  patterns  of  small  antennas  mounted  on  a  cylinder  (fuselage) 
can  be  found  using  the  previously  derived  solutions.  The  antennas 
considered  are  a  radial  ir.onopole  and  an  arbitarily  oriented  slot. 

All  have  dimensions  chat  are  small  in  terms  of  the  wavelength. 

The  radial  monopole  and  slot  are  considered  to  receive,  respectively, 
an  open  circuit  voltage  and  a  short  circuit  current,  so  that  they 
will  cause  negligible  distortion  to  the  field  which  would  exist 
there  with  the  antenna  absent.  Thus,  the  antenna  response  is  directly 
related  to  the  total  field  component  at  that  point. 

For  the  radial  monopole,  the  open  circuit  voltage  is  equal  to 
the  product  of  the  effective  height  of  the  monopole  times  the  component 
of  electric  intensity  parallel  to  its  axis  (E^).  For  the  slot  the 
short  circuit  current  is  equal  to  the  product  of  the  effective  height 
times  the  component  of  magnetic  intensity  parallel  to  its  axis  (cos  a 
H|  +  sin  a  Hj)  with  the  slot  orientation  defined  by  a  as  illustrated 
in  Fig.  2.  y 
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Fig.  2— Arbitrarily  oriented  slot  on  cylinder. 


It  is  shown  in  Ref.  1  that  the  effective  height  of  the  radial 
monopole  for  arbitrary  plane  wave  incidence  (e.,<j>.)  is  given  using 
Eqs.  (1)  and  (2)  by 

(3)  Monopole  =  Ce  E»p(p=b)  +  p  E>=b)] 


where  is  the  effective  length  of  the  monopole  without  the  cylinder 
present  and  b  is  the  distance  of  the  radial  monopole  from  the  center 
of  the  cylinder. 

The  effective  height  of  the  slot  can  be  written  similarly  as 

(4)  hslot  =  a  +  >CH^(p=a)  sin  *  +  H^2(?=a)c°si3}. 

Note  that  in  each  case  the  antenna  is  located  at  (;=0,  z=0)  and  the 
phase  is  referred  to  the  center  of  the  cylinder. 

The  transmitted  field,  when  the  antenna  elements  are  used  for 
transmission,  is[l] 
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Applying  the  above  solution  to  a  x/4  radial  monopole  it  was 
determined  that  the  series  converged  for  practical  purposes  (three 
significant  figures)  after  summing  (7a  sin  e-j  +  8)  terms.  The 
solutions  for  the  slot  antenna  converged  after  summing  (7.5a  sin  e.  + 
10)  terms. 

Since  these  expressions  for  the  distant  fields  radiated  are 
based  on  infinitesimal  elements  of  electric  and  magnetic  currents, 
the  far  fields  of  an  arbitrary  antenna  can  be  determined  by  integrating 
across  equivalent  aperture  currents.  This  is  'one  numerically  by 
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superposition  cf  discrete  elements  which  approximate  the  actual 
aperture  distribution.  This  technique  can  be  shown  to  be  valid 
provided  the  element  spacing  is  much  smaller  than  the  wavelength[22], 
that  is,  one  may  write 


(7)  jjtotal,  ^  ” 

1  1  n=l 


!n  hn^i,:i":n*zn* 


where  In  is  the  current  of  the  nth  element,  ( -n)  defines  its 
position,  and  hn  is  the  appropriate  effective  height.  For  example, 
some  of  our  slot  measurements  were  taken  using  open-ended  waveguide 
whose  aperture  distribution  can  be  adequately  approximated  by  six 
infinitesimal  sources.  Note  however  that  an  infinitesimal  antenna 
is  considered  in  each  case  except  where  it  is  stated  otherwise. 

These  solutions  (Eqs.  (1)  and  2))  are  not  just  limited  to  finding 
the  radiation  patterns  of  antennas  mounted  on  or  near  an  aircraft 
fuselage  as  used  previously.  They  are  applied  in  Chapters  III  and 
IV  to  determine  the  near  scattered  fields  by  aircraft  engines  which 
are  modelled  as  infinite  cylinders  in  Section  IV-B  and  finite 
cylinders  in  Section  III-C.  In  addition,  they  are  used  to  compute 
the  currents  on  an  infinite  cylinder  and  then  integrated  over  a 
finite  length  to  obtain  an  approximate  solution  for  the  far  field 
specular  scatter  of  a  finite  cylinder.  Consequently,  these  solutions 
are  used  extensively  throughout  this  report. 


C.  Wedge  Diffraction 

It  has  been  shown  in  the  last  section  that  a  cylindrically 
shaped  fuselage  can  be  analyzed  using  conventional  modal  solutions. 
However,  aircraft  shapes  are,  in  general,  quite  complex  involving 
many  complicated  structures.  Thus,  to  improve  the  radiation  pattern 
over  that  simply  obtained  using  the  modal  solutions,  one  must  take 
these  various  scattering  structures  into  account.  Based  on  past 
performance,  the  neometrical  theory  of  diffraction  has  proven  itself 
well  suitea  unis  type  of  analytical  study.  Not  only  does  this 
approach  fit  nicely  in  terms  of  the  ray  optics  format  but  it  also 
provides  a  means  for  analyzing  the  effect  of  three-dimensional 
structures.  Consequently,  a  brief  description  of  tl.e  wedge  dif¬ 
fraction  problem  is  presented  here.  The  diffraction  by  a  curved 
surface  is  treated  in  the  following  section.  Note  that  these  two 
basic  diffraction  problems  constitute  the  geometrical  theory  of 
diffraction  (GTD). 

An  asymptotic  solution  for  the  diffraction  from  a  conducting 
wedge  was  first  solved  by  Sonmerfeld[23].  Originally,  GTD[2]  as 
applied  to  diffraction  by  a  wedqe  was  based  on  plane  wave  diffraction 
coefficients;  however,  as  shown  in  Ref.  [25]  the  use  of  diffraction 
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of  cylindrical  waves  has  been  found  necessary  in  the  treatment  of 
antennas.  Consequently,  different  formulations  of  wedge  diffraction 
were  substituted  for  the  plane  wave  diffraction  coefficient  which  is 
the  basis  for  wedge  diffraction  theory.  Paulif26]  introduced  the 
Vr  function  as  a  prac;1cal  formulation  to  the  solution  for  a  finite- 
angle  conducting  wedge.  Recently,  however,  Hutchins  and 
Kouyoumjian[27,28]  presented  a  formulation  for  the  diffracted 
field  (Vg),  which  significantly  improves  the  accuracy  over  that  ob¬ 
tained  from  Pauli's  form. 

This  improved  diffraction  solution[27,28]  is  betcer  in  the 
transition  regions  (near  the  incident  and  reflected  shadow 
boundaries).  It  can  be  written  in  the  form 


(8)  VB(L,6,n)  =  I__(L,3-n)  +  I+_(L,c,n) 


where 

I±T7(L,6,n) 


,-j(kL+-/4) 

; — — - /a  cot  ( 

ja^T 


-±8 

“2n 


)  x 


ejkLa  f 


dt  +  [higner  order  terms] 


and  where  the  higher  order  terms  are  negligible  for  large  kL  and 
with  n  defined  from  the  wedge  angle  WA  =  (2-n )rr ,  also  a  =  1  + 
cos(8-2nTrN)  and  N  is  a  positive  or  negative  integer  or  zero,  which¬ 
ever  most  nearly  satisfies  the  equations 


2rorN-6  =  -tt  for  I 

-TT 

2nwN  -8=  +tt  for  I, 

+71 


The  variables  L  and  6  are  defined  later. 

The  three  dimensional  wedge  diffraction  problem  is  pictured 
in  Fig.  3.  A  source  whose  radiated  TT  field  is  given  by  IP(s)  is 
located  at  point  s'Cp'^'jZ1).  It  can  be  an  arbitrary  electric  or 
magnetic  source  causing  plane,  cylindrical,  conical,  or  spherical 
wave  incidence  on  the  wedge  tip.  The  diffracted  vector  field  at 
point  s(p,$,z)  can  be  written  in  terms  of  a  dyadic  diffraction 
coefficient.  Kouyoumjian  and  Pathak[29J  have  given  a  more 
rigorous  basis  for  the  GTD  formulation  and  have  shown  that  the 
diffracted  fields  may  be  written  compactly  if  they  are  in  terms  of 
a  ray-fixed  coordinate  system.  The  ray-fixed  coordinate  system 
is  centered  at  the  point  of  diffraction  Qg,  (or  points  of  dif¬ 
fraction  in  the  case  of  plane  wave  incidence).  is  a  unique 
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point  or  points  for  a  given  source  and  observation  point.  The 
incident  ray  diffracts  as  a  cone  of  rays  such  that  £0  =  (see 
Fig.  3). 

The  relationships  between  the^orthogonal  unit  vectors  associated 

with  these  coordinates  (s' , 2  .:)  are  given  bv 

0  0 

I  =  -  s' 

I  =£  'x 
o 

a  a  <«*» 
s  =  5o  x 

where  I  is  the  incident  direction  unit  vector,  and  s  is  the  dif¬ 
fraction  direction  uni*  vector.  The  diffracted  fielo  is  now  given 
by 


Hs)  =  e1(qe)  •  De(s,I)  A(s)  e“3ks. 

For  our  purpose,  it  is  more  convenient  to  write  the  diffracted  field 
in  terms  of  the  Vg  function  in  Ee.  (1)  as,* 


E”(s)l  [  ~vb  °  i  f 


E^(s) 


iU-  a(s)  e-Jks 


sin  ft 


o  -v+J  le2(qe) 


where 


Vg  =  VB(L,8",n)  +  VB(L,8+,n). 

The  minus  sign  (Vg)  applies  for  the  E-field  component  parallel  to 
the  edge  with  boundary  condition 


^ wedge ^  "  °* 

The  plus  sign  (Vg)  applies  for  the  "E-field  vector  perpendicular  to 
the  edge  with  boundary  conditions 


*If  a  fixed  coordinate  system  is  used  Eq.  (9)  takes  the  form  of  a 
3x3  matrix. 
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The  annular  relations  are  expressed  by 


i~s  =:  h' 

where  the  minus  sign  (b~)  is  associated  with  the  incident  field  and 
the  plus  sign  (5+)  with  the  reflected  field.  The  quantity  A(s)  is 
a  ray  divergence  factor  given  by[29] 

1  plane,  cylindrical  (s=p), 

and  conical  wave  incidence 


A(s)  =< 


s' 

hW- 


+s) 

ana  L  is  given  by[29] 


spherical  wave  incidence 


f 

s  sin28„ 
o 


L 


J  p'p 

o+o ' 


2 

s's  sin  Bq 
s+s‘ 


plane  wave  incidence 


cylindrical  wave  incidence 


conical  and  spherical  wave 
incidence. 


For  the  two-dimensional  wedge  problem,  illustrated  in  Fig.  4, 
where  there  is  cylindrical  wave  incidence  with  8  =  90°,  Eq.  (9) 
reduces  to  give 


E^(p,4>) 


e;,(p  •,*•)' 
E> 


Jkp ' 


N  p+p' 


In  the  far  field  (p>>p‘)  this  becomes 


V 


1 


OBSERVATION  POINT 


Putting  t';is  in  ray  form  and  factoring  out 


00) 


-V-  o' 

Rj(o’)' 

-  0  -VB. 

we  have 


The  ray  form  used  here  is  given  by 

P-o‘kp 

E(p,$)  =  RU)  e— —  • 

4> 

Thus,  R’U)  is  related  to  the  far  field  pattern  function. 
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For  the  three-dimensional  wedge  problem  where  there  is  spherical 
wave  incidence.  Eq.  (9)  reduces  to  give 


E||(s.60>;) 

r 

O 

1  CO 
> 

1 

E?(S,S0,:)_ 

-VB- 

L 

In  the  far  field  (s>>s'),  we  have 


s_ _ 

s+s‘ 


2 

tfks's  sin  B. 


s+s 


')* 


ks 


’e?(s,b0>S)' 

'-V-  o  ■ 

r 

-eJ(s,20,:)_ 

-°  -  VB. 

2 

jks'sin  £  -jks 
s'e  0  ^-r— 


e-jks 

Putting  this  in  ray  form  and  factoring  out  — -  ,  it  is  seen  that 


Hi) 


i2<v>> 

•£<V*> 


-jks'  cos^P 
e  iJ 


It  is  interesting  to  note  that  in  the  principal  plane  (e0  =  90°)  the 
ray  form  of  the  three-dimensional  case  takes  on  the  same  form  as 
the  two-dimensional  problem. 

The  total  ray  value  at  the  observation  point  s  is  given  by  the 
sum  of  the  geometrical  optics  terms  and  the  diffracted  terms 


RT(s)  =  RG*°*(s)  +  R^s) 


where 


R6*0* 


(s) 


f-i 

1  R  (s)  +  R  (s)  incident  and  reflected  region  I 
-  R°(s)  incident  region  II 

0  shadowed  region  III 


and  R  (s)  may  be  determined  from  the  image  of  the  source  term  using 
basic  geometrical  optics  techniques.  These  three  regions  are  il¬ 
lustrated  in  Fig.  4  for  a  two-dimensional  wedge  diffraction  problem. 


D.  Curved  Surface  Diffraction 


When  an  incident  ray  strikes  a  smooth,  curved  perfectly  con¬ 
ducting  surface  at  grazing  incidence,  i.e.,  at  the  shadow  boundary, 
a  part  of  its  energy  is  diffracted  into  the  shadow  region.  To 
describe  this  phenomenon  Keller[30l  introduced  a  doss  of  curved 
surface  diffracted  rays.  These  ray  paths  include  the  points  Qi  and 
Q2  which  form  a  curve  on  the  diffracting  surface  as  illustrated  in 
Fig.  5.  However,  the  actual  concept  of  creeping  waves  was  introduced 
by  Franz  and  Deppeman.L31 ,321  The  basic  concept  as  presented  in 
the  following  discussion  is  basically  taken  fr cm  "Asymptotic  High- 
frequency  Methods"  by  Kouyoumjian.I33~ 


SHADOW 

B0UND4RY 


Fig.  5— Diffraction  by  a  smooth  curved  surface. 


The  diffraction  by  a  smooth  curved  surface  is  shown  in  Fig.  5 
in  which  0  is  the  source  point  and  P  is  the  observation  point  in 
the  shadow  region.  Applying  Fermat's  principle,  the  line  OQi Q2P 
is  the  shortest  distance  between  0  and  P  which  does  not  penetrate 
the  surface.  In  detail,  a  ray  incident  on  the  snadow  boundary  at 
Ql  divides;  one  part  of  the  incident  energy  continues  straight  on 
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as  predicted  by  geometrical  optics,  a  second  part  follows  the 
surface  s  into  the  shadow  region  as  a  surface  ray  shedding  dif¬ 
fracted  rays  tangentially  as  it  propagates.  With  t,  n  and  B  being 
the  unit  vectors  in  the  direction  of  incidence,  normal  to  the 
surface  s  and  binormal  to  the  surface  (B  r  t  x  n),  respectively,  the 
incident  field  P_{0i)  may  be  resolved  into  its  normal  and  tangential 
components  (n  •  EMQl)  and  B  •  EMO]}).  It  is  assumed  that  these 
two  components  induce  surface  ray  fields  which  propagate  inde¬ 
pendently  of  each  other  along  the  geodesic  arc  between  Q]  and  Q2. 

From  Reference  1331  the  surface  ray  field  AeO;  at  C,  is  related  to 
the  incident  fields  at  Q-j  by 

(12)  A(0,)  e  1  =  Ds(n])  b,  .  E1(f)1> 

where  Ds(Qi)  is  the  scalar  diffraction  coefficient  for  a  soft  surface. 
The  amplitude  of  the  surface  ray  is  assumed  to  be  governed  by  the 
conservation  of  energy  between  a  pair  of  adjacent  surface  rays. 

Hence,  the  amplitude  behavior  of  the  fields  is  given  as 

K  f-.v(t')dt'] 

(13)  A(Q2)  =  A(Q-j)  e 

where 

dn-j  and  dn2  =  the  separation  between  a  pair  of  rays  at 
Q-j  and  Q2,  respectively. 

a(t)  =  the  attenuation  constant  which  is  a  function 
of  t,  the  coordinate  along  the  surface  ray, 
because  it  depends  on  the  local  radius  of 
curvature  and  its  derivatives. 

The  attenuation  constant  a(t)  is  introduced  due  to  the  tangential 
shedding  of  rays  as  the  surface  ray  propagates.  It  is  seen  from 
Fig.  5  that  Q2  is  a  caustic  of  tie  diffracted  field  and  the  second 
caustic  is  located  at  a  distance  p  from  Q2.  Thus,  the  tangential 
component  of  the  diffracted  field  which  radiates  from  Q2  towards  P 
can  be  found,  as  in  the  previous  edge  diffraction  case,  with  one  of 
the  caustics  used  as  a  reference  point  and  is  given  by 


(14) 


eV) 


ds(q2)  a(q2) 


04>(Qo) 


/ 


s(p+s) 


,-jks 
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From  Eqs.  (12),  (13)  and  (14)  there  results 


(15) 


b2  ■  fw  =  [s,  •  e’iq,)]  Ds(0t)  os(n2)J  §  Lfcj 
r  i°2  I 

-3  k(t+s)  +  ;  •:* ( t * )  dt ' 

L  Q, 


It  is  found  that  b-j  •  E1(Qi)  excites  an  infinity  of  surface 
ray  inodes  each  with  its  own  diffraction  coefficient  and  attenuation 
constant.  Thus,  the  expression  in  Eq.  (15)  is  replaced  by 


(16) 


b9  •  EV) 


dr>i 


S (n+S ) 


s-j  [k(t+s  jj 


L  m 


7  D  (Q,)  D  (Q„)  e  -1 
£  smvvl '  smw2' 


r  2 

-i  t  (t* ) 

>n  mv 


dt' 


Equation  (16)  relates  the  diffracted  field  at  P  to  the  incident  field 
at  Qi  for  the  soft  surface  boundary  condition. 

An  expression  similar  to  Eq.  (16)  is  also  obtained  for  the 
normal  component  of  the  incident  field;  in  this  case,  the  scalar 
diffraction  coefficients  and  attenuation  constants  for  the  hard 
surface  replace  those  of  the  soft  surface.  Therefore,  the  vector  dif 
fracted  field  at  P  can  be  written  in  terms  of  the  electromagnetic 
field  incident  at  O-j  as 


(17)  E(P)  =  [n^  v(l  ,2)  +  b2b1  u(l  ,2  j]  •  E^l  )Jj 
in  which  v(l,2),  u(l,2)  are  equal  to  0 


-jks 


(18) 


dn-i  «  "J 

_L  e  Jkt  7  D  (1 )  D  (2)  e 
ario  £  m  m  ' 

c  m 


am(t')dt' 


in 


with  the  subscripts  h,  s,  respectively,  added  to  Dm  and  am.  Note 
that  Q-j  and  Q2  have  been  replaced  by  1  and  2  for  the  sake  of  brevity. 
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Finding  dr-j,  d^2>  and  p  is  simply  a  matter  of  differential  geometry 
involving  the  rays  and  surface;  this  is  discussed  at  length  in  Levy 
and  Keller[303.  The  generalized  diffraction  coefficient  and 
attenuation  constant  can  be  found  in  Reference  [343. 

The  diffraction  thus  far  discussed  is  applied  to  the  open 
curved  surface.  For  a  closed  surface,  each  surface  ray  mode 
produced  at  Q]  encircles  the  surface  an  infinite  number  of  times. 
The  length  of  the  surface  ray  path  for  the  n-th  encirclement  is 
t  +  nT  where  T  is  the  circumference  of  the  closed  surface.  These 
multiple-encircling  rays  may  be  summed  to  contribute 

rT  3 

f  "jkT-|  am(f)df 
\  1  ' 0 


to  the  denominator  of  the  diffracted  field.  It  is  interesting  to 
note  that  there  must  be  another  pair  of  diffraction  points  Q3  and 
Q4  for  the  closed  surface  as  shown  in  Fig.  6.  Therefore,  the  field 
at  any  point  P  in  the  shadow  region  is  the  sum  of  these  two  dif¬ 
fracted  fields  frcm  Q]-Q2  and  Q3-Q4.  The  total  field  at  any  point 
in  the  illuminated  region  is,  by  the  superposition  principle,  the 
sum  of  incident,  reflected  and  diffracted  fields.  A  detailed 
discussion  of  this  subject  can  be  found  in  Reference  [343- 

An  important  special  case  of  this  GTD  solution  is  the  one  in 
which  the  antenna  is  mounted  directly  on  the  curved  surface.  This 
problem  has  been  studied  by  Kouyoumjian[353  in  which  he  analyzed 
by  asymptotic  methods  the  far  field  patterns  of  various  antennas 
mounted  on  an  infinitely  long  perfectly  conducting  circular 
cylinder.  In  this  solution  the  reciprocity  theorem  is  employed  so 
that  the  radiation  pattern  is  derived  from  the  plane  wave  scattering 
by  a  cylinder.  The  physical  optics  approximation  for  the  surface 
current  is  used  to  describe  the  field  in  the  region  in  which  the 
source  is  directly  illuminated  by  the  incident  plane  wave  (the  lit 
region).  The  Fock  approximation[363  is  used  for  the  penumbra 
(transition)  region;  whereas,  in  the  deep  shadow  region  of  the 
cylinder  the  geometrical  theory  of  diffraction  solution  is  applied 
to  give  the  far  field  expression.  It  is  noted  that  a  launch  coeffi- 
fient,  which  relates  the  creeping  wave  (GTD  field)  to  the  actual 
surface  field  or  current,  is  introduced  so  that  the  GTD  solution 
can  be  properly  employed.  This  launch  factor  is  deduced  from  the 
exact  solution  for  the  surface  current  on  the  circular  cylinder;  and, 
in  its  asymptotic  form,  involves  the  well-known  Airy  function. 

The  far-zone  principal  plane*  pattern  for  an  infinitesimal 
slot  mounted  on  an  infinitely  long  perfectly  conducting  circular 


*The  principal  plane  being  a  plane  cutting  across  the  cross-section 
of  the  cylinder. 
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[5^1*] 


J; 


S3 


LOCATION  OF  THE  SOURCE 
AND  DIFFRACTION  POINTS  A ,8 
A,  B 


REGION  I  :  LIT  REGION 
REGION  Ti:  TRANSITION  REGION 
REGION  EX'  DEEP  SHADOW  REGION 


Fig.  7— Antenna  mounted  on  a  circular  cylinder:  GTD  solution 
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where 


•1  "  -o 


'z  =  T  "  !6o 
M  =  (  fi)1/3 


r  -  (  kaxl/3  , 

C1  ( T*  *1 


r  _  /  ka  xl/3 

:r  v2 


-/2<e<3-/2 


The  functions  g*(ci)  and  g*{£2)  are  the  conjugate  of  g(r])  and  q(ro) 
respectively.  The  function  g(c)  is  the  so  called  Fock’s  function' 
for  hard  boundary  condition  (au/'*n  =  0)  and  is  given  as 


9(0  = 


where 


Jr,  ^  dT  ’ 


Wo(t) 


tt  Jr 


eTt‘  ™  dt 


and  r-]  and  r2  are  contours  in  the  complex  plane  as  shown  in  Fig.  8. 
The  magnitude  and  phase  of  the  Fock  function  [ofxll  are  qiven  in 
Reference  [36].  ‘ 

3)  Deep  Shadow  Region:  Surface  Ray  Description 

_l/3  ■;  v 

(2>)  RyCo.e)  =  (-2,)  (jji)  e  12  • 

I  <0*0*  S'  e'JkS' 

m=  o 


xh/  »-|2  ft •/  ~  \  "Cam(a)^s2  ~jks2 
L(a)3  Ai(-q  )  e  e  1  } 


•r  X 


*  J 


Fig.  8— Contours  for  complex  plane  integrations. 


where 

s-j  =  a^ 

s2  -  a ^ 

and  and  xpo  were  defined  earlier  and  %  is  the  root  of  Ai1  (-%), 
that  is,  Ai'[-qm)  =  C.  The  diffraction  coefficient  [Dl](a)l2  and 
the  surface  attenuation  constant  ol](a)  are  given  in  Taole  I.  The 
roots  qm  and  AI(-qm)  are  given  in  Table  II.  It  is  noted  that  the 
surface  ray  expres  ion  in  Eq.  (21)  can  be  rewritten  in  the  following 
way 


(22)  R,(p,e)  =  (-2n)  (  fj 


ka, -1/3  T? 


I  WlWf  Ai(-qje 


-Y(a)si 


m=o 


+  rn^a'i2 


D>m<a>r  Ai(-qm)e 


_  -y(a}s; 


where  y(a)  -a  ^(a)  +  jk;,  propagation  constant  of  the  surface  ray. 
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TABLE  II 

AIRY  FUNCTIONS  AND  ROOTS 


BOUNDARY  CONDITION 

MODE 

SO 

FT 

HARD 

m 

A.(-q  ) 
iv  Hnr 

°m 

0 

2.33811 

0.70121 

1.01879 

0.53566 

1 

4.08795 

-0.83011 

3.24820 

-0.41902 

2 

5.52056 

0.86520 

4.82010 

0.38041 

3 

6.78671 

-0.91085 

6.16331 

-0.35791 

4 

7.94413 

0.94734 

7.37218 

0.34236 

The  calculated  radiation  patterns  using  Eqs.  (19)  to  (22)  in  the 
xz-plane  (principal  plane)  of  Fig.  7  for  various  cylinders  have  been 
shown  to  agree  very  favorably  with  those  obtained  from  the  modal 
solution,  especially  for  ka  >_  3.0.  These  formulations  also  predict 
the  radiation  pattern  fairly  well  for  small  cylinders  with  radii  down 
to  the  range  from  1/3  to  1/4  wavelength. [37] 

The  results  for  the  circular  cylinder  are  extended  to  cylinders 
with  general  cross-sectional  shape  in  which  the  curvature  varies 
along  the  ray  path.  This  variable  curvature  has  a  strong  effect  on 
the  attenuation  constant  which  in  turn  affects  the  energy  propagation 
quite  significantly.  Thus,  the  representations  for  fields  in  the 
lit,  transition  i  id  deep-shadow  region  are  modified  to  include  the 
variable  curvature  effect  and  again  are  given  by 

1)  Lit  Region: 

jkr  cos(e-en) 

(23)  Ry(e)  =2e  p 


2)  Transition  Region:  p 


(24) 


kds 


-jf 


kds 


Rw(e)  =  {g*(cfl)  e  A  •  F(A)  +  g*UR)e  ;B 


F(B)} 
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3)  Deep-Shadow  Region: 


T  k‘a  ~1/31 

(25)  Ry(  )  =  -2-j|(  y-2-) 


eJ  3  (  ^  )V2, 


-Ji  >m(s)  ds 


•  UD>>  D>  «<-%,> e  'A 

m=o  I 

-3 !  ds 

+  Dr^(P)  D^(B)  Ai  (— qV )  e  'B 
rn  m  m 


F(B) 


F  (A) 


where  A  and  B  are  the  points  that  the  surface  rays  diffract  toward 
the  direction  of  the  observation  point,  «;  (rs,6p)  defines  the  source 
locations;  and  p  is  the  radius  of  curvature  of  the  curved  surface; 


and 


f 

-'b 


V  X1/3 

Itf)  ds 


The  increment  of  the  arc  length  along  the  ray  trajectory  is  ds; 

F(A),  F(B)  are  the  phase  correction  factors  at  the  tangent  points 
A  and  B;  g*U/\),  g*( €b)  are  tkie  conjugates  of  gUa)  and  9Ub)>  the 
Fock  function  as  defined  previously;  and  Dh(P),  D'j’(A),  and  D^(B) 
are  the  diffraction  coefficients  at  the  source  location,  and  points 
A  and  B,  respectively.  The  propagation  constant  is  y ^ ( s )  =  o£(s)  +  jk, 
where  a ”  is  the  ray  trajectory  attenuation  factor  for  the  hara  boundary 
case  ana  k  is  2n/\.  The  general  expressions  used  for  the  diffraction 
and  attenuation  constants  are  given  in  Table  I.  Note  that  the  phase 
reference  point  for  these  formulations  is  at  the  origin.  Some  calcu¬ 
lated  radiation  patterns  in  xz-plane  (the  elevation  plane)  employing 
these  formulations  for  an  elliptical  cylinder  are  also  given  in 
Reference  [37],  These  results  compare  fairly  well  with  the  exauc 
solution  obtained  by  Sinclair[l]  as  shown  in  Fig.  9. 
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NORMALIZED  MAGNETIC  FIELD 


\ 


Fig.  9--Principal  plane  patterns  with  a  circumferential 
slot  mounted  on  an  elliptical  cylinder. 


The  same  technique  is  also  employed  to  derive  the  generalized 
expressions  for  a  slot  antenna  with  soft  boundary  condition  and  in¬ 
finitesimal  monopole  antenna.  The  equations  again  contain  the  Fock 
functions.  A  sunmary  of  the  far  field  expressions  for  the  circumfer¬ 
ential  slot,  axial  slot  and  infinitesimal  monopole  are  presented  in 
Table  III[35]. 
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CHAPTER  III 

NEAR  FIELD  SCATTERING  BY  FINITE  STRUCTURES 


A.  Introduction 


Two  relatively  new  near  field  scattering  problems  are  considered 
in  this  chapter  using  various  high  frequency  approximations.  In 
these  problems  the  antenna  is  in  the  near  field  of  a  finite  three- 
dimensional  scattering  body;  whereas,  the  observation  point  is  in 
the  far  field.  However,  it  is  assumed  that  any  point  on  the 
scattering  body  s  in  the  far  field  of  the  source.  Thus,  if  a  1  rge 
antenna  illuminates  one  of  these  scattering  bodies,  one  must  in¬ 
tegrate  over  the  equivalent  aperture  currents  using  the  equivalence 
principle  to  obtain  the  currents  and  the  radiation  integrals  to 
compute  the  far  field  pattern.  However,  as  shown  in  Ref.  [22],  this 
integration  can  be  accomplished  numerically  using  a  finite  number 
of  discrete  infinitesimal  sources.  This  point  was  discussed  in 
Section  II-B. 

The  scattering  body  is  defined  in  terms  of  its  location  and 
dimensions  in  a  general  coordinate  system,  and  the  antenna  is  defined 
by  its  location  and  far  field  pattern.  Note  that  some  assumptions 
are  made  in  terms  of  the  geometry  so  that  these  r  ~ults  ultimately 
find  application  in  terms  of  our  approximate  on-aircraft  antenna 
problem. 

These  analyses  are  applications  of  the  fundamental  problems 
considered  in  the  previous  chapter.  Similar  types  of  studies  have 
been  made  previously  such  as  the  near  field  scattering  by  various 
geometric  bodies  considered  by  Lentz[38].  Near  field  scattering 
effects  must  be  included  in  an  analysis  of  radomes  such  as  the  plane 
wave  spectrum  analysis  done  by  Wu[39].  Thus,  near  field  scattering 
problems  are  of  recent  interest  and  can  have  application  beyond  the 
on-aircraft  antenna  problem  considered  here. 

B.  Near  Field  Scattering  by  a  Finite  Plate 

The  near  field  scattering  by  a  finite  flat  plate  is  a  relatively 
new  topic  at  higher  frequencies  where  the  plate  is  large  in  terms 
of  the  wavelength.  The  solution  presented  here  is  a  practical  appli¬ 
cation  of  the  three-dimensional  we dge  diffraction  theory  given  earlier. 
The  flat  plate  geometry  is  illustrated  in  Fig.  10.  Tne  source  is 
defined  by  its  location  and  far-field  pattern.  The  far-field  pattern 
of  the  source  is  appropriate  in  that  the  plate  is  located  at  least 
2D 2/x  away  from  the  source  where  D  is  the  maximum  dimension  of  the 
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I  l«l.  '(!  Him  t  I told  f  lat  pi  ate  geometry. 

.•inn  n  iim  finite  M.it  plate*  is  simply  specified  by  the  location  of 

I  I  ■<  II  I  III  lllil  S  . 

l(  i*.  hniwn  fli.i  1  fin  n  given  scatter  direction  there  is  only 
mm  i >ii int  tt 1 1*11*1  .in  tut  tn 1 1 el v  long  straight  edge  at  which  the  dif- 
hattis!  r  1 1* I .(  v  .in  emanate  tor  a  near  zone  source.  Thus,  this  point 
must  In-  tmiiiil  tr  it  ovists  tor  each  of  the  n  edqes  that  describe 
t  lir  flat  p  I  rt  t  e . 

lot  ns  new  , ensi.ter  rt  general  nth  edge  using  the  geometry  il- 
lnst*ate*t  i"  tig  tl.  the  e\e**e*sion$  necessary  to  f’nd  the  cir- 
“rt.tien  is»iiit  e  given  tv! On : 
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Fig.  11— Geometry  for  diffraction  point  along  the  mth  edge. 


It  is  known  that  the  angle  of  incidence  (3^)  is  equal  to  the  angle 
of  diffraction  ( 80) ,  which  in  turn  is  defined  by  the  given  scatter 
direction  which  is  given  by 

d  =  sin  e$  cos  <^sx  +  sin  eg  sin  <*>sy  +  cos  esz  . 

Since  these  angles  are  equal  (s'  =  b  )  the  cosines  of  these  angles 
must  be  equal  which  implies 


d  •  e  = 
m 


sin  es cos  :s(VrV  ,  sin  As  sin  :s(Vi^> 


COS  6  (z  , -Z  ) 
sv  m+I  m' 

ft 


Thus,  one  finds  that 


x  »  t(yry .  \  (  ,  j  .towi  „ 

m  H  xs  lVl  V  ym  M  y: 


(Vr-vm>  +  An  *  t(-ZnMl  Zm>  -  0  ( WO  *  A, 


s  >vrv  -  *■< 


‘  -  jf  CAv-(xn,-xs)(xn,+rxin)-(y|n-y5)(y|n+r  »J-(v2s)(W 


where 


H  N'Vr*/  +  <VrV2  +  (zm*rzm>2 


sin  «s  cos  ♦sfVl-V^i^s  sin  *s(Vl-».,tcos  es!Wz» 


.1-0  V  /  t(ym,,-ym)  v 


An  iteration  technique  can  now  be  applied  using  the  above 
equations  by  placing  the  present  (or  assumed)  value  of  t  in  the 
above  equation  for  A  and  then  proceed  to  solve  for  a  new  value  for 
t  using  Eq.  (28).  It  can  be  shown  analytically  that  this  iteration 
approach  will  always  converge  to  the  desired  value  of  t,  which 
locates  the  point  along  the  edge.  However,  iteration  techniques 
often  are  not  efficient  in  numerical  computations  as  was  found 
in  this  case.  It  was  found  that  in  certain  regions  too  many 
iterations  were  required  which  resulted  in  a  loss  of  numerical 
accuracy.  Consequently,  a  search  approach  has  been  developed  in 
which  the  values  of  esdp  *  em  are  determined  at  the  end  points  of 
an  edge  of  the„plate  anci  at  its  midpoint.  These  values  are  then 
compared  with  d  *  em  for  which  one  can  easily  determine  whether  the 
diffraction  point  falls  within  the  limits  of  the  finite  edge  and 
if  so  which  side  of  the  midpoint.  These  tests  are  all  based  on  the 
dot  product  comparisons.  For  example,  if  the  value  of  d  •  em  does 
r.ot  fall  between  the  values  of  escjp  •  em  for  the  two  end  points, 
then  a  diffraction  does  not  occur  from  the  mth  edge  for  that  scatter 
direction.  If  the  scatter  direction  dot  product  falls  between  an  end 
point  and  the  mid  point,  ont  then  determines  e$dp  *  em  for  the  midpoint 
of  the  new  region  within  which  a  diffraction  occurs.  From  these  values 
one  finds  a  smaller  region  within  which  a  diffraction  must  occur.  This 
process  continues  until  the  diffraction  point  is  found  v/i thin  certain 
minimum  limits. 

Once  the  diffraction  point  is  located,  one  must  find  the  dif¬ 
fracted  field  value  from  the  mth  edge.  The  far  field  pattern  of  the 
source  can  in  general  be  written  as 

— j  ks 1  — jks* 

(29)  E$(0,4>)  =  [0  F(e,>f)  +  $G(e,<j>)U  —  =  $(0,$)  — 

where  s'  is  the  range  from  the  source  to  the  field  point.  Using  the 
geometry  illustrated  in  Fig.  12a  and  applying  the  results  presented 
in  Section  Il-C  one  finds  that 


where 

R  u  ~  ^  ( 0  j  >  4*  -j )  * 

R1  =  R(0-j  »4>.j)  * 
2 

kpp=  ks'  sin  $Q 
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Fig.  12a— Edge  diffraction  geometry. 


Y  =  xdp  sin  es  cos  +  ydp  sin  e$  sin  *s  +  Z(Jp  cos  es 

VB  =  VB'kpp>  *  VB^pp’  *+V2)' 

The  incident  field  direction  (e.,4>.)  from  the  source  to  the  point  of 
diffraction  is  defined  by 
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Note  that  these  angles  as  well  as  all  others  are  defined  in  terms  of 
an  arctangent  in  that  for  numerical  results  the  Fortran  arctangent 
routines  define  the  angle  between  -180°  and  180°.  It  is  easily  shown 
that 

sin  B0  =  Id  x  ej 

where 


0°  <  60 <  180°. 

In  order  to  find  the  diffraction  angles  it  is  necessary  to  define  a 
coordinate  system  at  the  diffraction  point  as  illustrated  in  Fig.  12b. 


Fig.  12b— Edge 


“fint  of  diffractir“ 


A 


>IQI« 


ill] 


The  three  orthogonal  unit  vectors  are  e^  thermit  vector  along  the 
edge,  n  the  normal  to  the  surface,  and  ejj,  =  n  x  The  incident 
vector  direction  can  now  be  written  in  terms  of  this  coordinate 
system  as 


-I  =  e1  sin  co°  a  +  n  sin  e  sin  e  i-  e„  cose„. 
moo  o  o  m  o 

Using  the  above  expression  one  finds  that 


The  scatter  direction  unit  vector  can,  also,  be  decomposed 
in  terms  of  this  coordinate  system  by 

A  A  A  A 

d  =  e^  sin  e  cos  4  +  n  sin  e  sin  v  +  em  cos  e 
which  results  in  giving 


The  vector  directions  of  the  diffracted  field  are  defined  by  the 
following  expressions: 

A  ^  ^ 

4>  =  -e'  sin  ®  +  n  cos  4> 
o  m  o  c 

4>=  -e'  sin  <t>  +  n  cos  4> 
m 

<*•  A  A  A  A  A 

e;  =  *o  x  1  and  6q  =  «i>  x  d. 

Once  these  terms  are  determined  the  total  diffracted  field  from  a 
general  mth  edge  is  given  {Eq.  (30))  by 

(31)  R^Ces,4.s)  =  R^s0  +  R^$  . 
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Using  the  superposition  principle  the  total  singly  diffracted  field 
by  the  n  edges  of  the  flat  plate  is  giv°n  using  Eq. (31 )  by 

(52)  S’fvV  -  X  **  (V4S>  ■ 

The  reflected  field  from  the  flat  plate  is  considered  next.  The 
first  step  in  this  calculation  is  to  find  the  location  of  the  image 
source,  whicn  is  uniquely  determined  once  the  plane  of  the  flat  plate 
is  defined  relative  to  the  source  location.  In  fact,  the  image  is 
located  along  a  line  which  is  orthogonal  to  the  plate  and  positioned 
an  equal  distance  on  the  opposite  side  of  the  plate.  This  location 
can  be  found  analytically  using  the  geometry  illustrated  in  Fig.  13a. 


SOURCE 

LOCATION 


Fig.  13a.  Geometry  for  location  of  image  source. 
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The  image  position  vector  (R.)  can  be  determined  using  the  .^ollrwing 
expressions: 


(33) 


R.  =  - 1  n  =  2{"r,  +  t-.  t-i  +  t^to  1 
n  n  t  ]  ll  22J 

R-j  =  (xrx$)x  +  (y-,  -ys  )y  +  (z-,-zs)z 


n  =  nxx 


+  n  y  + 
T 


V 


=  V  +  V  *  ^z2 


l2  =  l2xx  +  V  +  ^z2 


which  result  in  the  following  matrix  equation 


(-  > 


Note  that  n,  t-| ,  and  to  are,  respectively,  the  unit  normal  and  tangent 
vectors  associated  with  the  plane  of  the  flat  plate.  The  above  matrix 
can  now  be  inverted  giving  the  values  of  ciq,  a],  and  which  in  turn 
defires  fL.  Knowing  the  vectcr  ( Rn) ,  the  image  vector  location  is 
simpl  •/  defined  using  Eq.  (33)  by 


where 

A  A 

Rs  =  xsx  +  V  +  zsz  • 

With  the  image  position  known,  one  needs  to  determine  if  the 
reflected  field  contributes  to  the  total  scattered  field  using  the 
geometrical  optics  approach.  If  the  reflected  field  is  a  con¬ 
tributor,  the  ray  from  the  image  source  in  the  scatter  direction 
must  pass  through  the  finite  plate  limits.  Thus,  one  must  find 
the  location  of  the  intersection  point  of  this  ray  and  the  plane 
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containing  the  flat  plate.  This  is  accomplished  using  the  geometry 
illustrated  in  Fig.  13b.  The  relationship  necessary  to  find  this 
point  is 

Yl*l  +  Y2*2  “  V  =  *i  "  ^1 
which  results  in 


\x 

*2x 

-sin  cos 

;s  ) 

1  ' 

/xi  -  X1  \ 

t2y 

-sin  es  sin 

A 

vs 

yi  -  yl 

t2z 

-cos  &s 

) 

W 

\zi  "  Z1  1 

Again  the  above  matrix  can  be  inverted  to  solve  for  the  point  of 
intersection. 


SOURCE 


Fig.  13b— Intersection  of  image  ray  with  flat  plate. 
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One  must  now  decide  if  this  point  falls  within  the  finite  limits 
of  the  flat  plate,  which  can  he  accomplished  many  ways.  The  scheme 
used  here  is  simply  to  determine  if  the  angle  (-mi)  between  one  of 
the  edges  of  a  general  corner  and  a  line  going  from  the  general 
corner  to  the  intersection  point  is  greater  than  the  angle  (  m)  between 
the  edges  making  up  the  corner.  If  is  greater  than  for  any  of 
the  corners,  then  reflections  do  not  occur.  Examples  of  this  technique 
are  shown  in  Fig.  14.  Note  that  reflections  occur  for  the  geometry 
of  Fig.  14a  since  rm-j  ■-  em  for  all  m;  whereas,  reflections  .)  not 
occur  in  Fig.  14b  since  emi-  >  at  the  mth  corner. 

Once  it  is  determined  that  reflections  do  occur,  it  is  necessary 
to  find  the  value  of  the  reflected  field.  It  is  known  that  the 
reflected  field  from  the  image  source  can  be  written  as 


(34) 


:rtr(es. 


using  the  geometry  illustrated  in  Fig.  15.  The  ray  form  of  the 
reflected  field  is  given  simply  by 


jk[  x.  sin  o$  cos  $  +  yi  sin  *s  sin  >  +  zi  cos -i  J 
e 


with  er  and  4>r  being  related  to  the  image  source  coordinate  system. 
The  above  quantities  can  be  found  from  the  boundary  conditions  that 
must  be  satisfied  on  the  flat  plate  and  which  are  given  by 


(36)  n  •  Er  (on  plate)  =  n  •  Es  (on  plate) 

t-j  •  Er  (on  plate)  =  -t-j  •  Es  (on  plate) 

•  E  (on  plate)  =  •  E  (on  plate). 

->r  +s 

Note  that  E  and  E  are,  respectively,  the  image  and  source  field 
values  at  the  point  of  intersection  on  the  flat  plate.  The  geometry 
used  to  define  these  various  terms  is  shown  in  Fig.  13b  with 


Rr  =  Rt  -  Rs  =  (xt  -  xs)J  *  (yt-ys)y  +  <vzs)z 
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Fig.  14a— Reflections  occur. 


(b) 


Fig.  1 4d — Reflections  do  not  occur. 
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Fig.  15— Reflected  field  geometry, 
and 

R  =  sin  e.  cos  $.  x  +  sin  e.  sin  y  +  cos  e-z 
r  i  i  i  1  1 

which  results  in  the  incident  field  direction  defined  by 


Thus,  at  the  point  of  intersection  one  finds  that 


-jkrr 


Es  (on  plate)  =[ebF(ei-, t.)  +  ::>G(vi-,;i)] 


-~s 


(37) 


with 


Pr 


Er  (on  plate)  =[erFr(er, :r)  +c'G' (?r,:r)Z 


r„r, 


-•krf 


ri 


r  =  r-.  Substituting  the  above  expressions  in;o  Eq.  (36)  one 
fir.dssthat1 

n-Es  =  [Fes*x  +  G$s*xZn  +  LFeS*y  +  G;s.yj  n  -  [Fes.z  + 

X  j 

G^s‘z]  nz 

n .r  =  CFr0r-x  +  Gr7-x X  +  CFrer*-y  +  Grir-J  ]n 

x  y 

y*A  a  a  a  A 

+  ZF  e  •  z  +  G  j  *z]  n 


with 


es  =  x  cos  e.  cos  +  y  cos  e.  sin  <>.  -  z  sin 

a  q  A  A  A 

<T  =  -x  sin  $  .  +  y  cos 

er  =  x  cos  es  cos  +  y  cos  es  sin  4>s  -  z  sin  es 


=  -x  sin  4>s  +  y  cos  $s  . 


A  A  -*S 

One^can.^alsQ,  write  similar  equations  for  ((x-)t]+a2t2)*E  and 
(a,t-+a9t9)-t“.  The  resulting  equations  can  be  suirmarized  by 


*ri  “2k2 

S  •  Er  -  F1  C„  +  6‘  C 


12 


n  •  E  =  F  C-|2  +  G  C-j^ 

(a-|  t*j  +  a^2 )  •  Er  =  Fr  C2^  +  Gr  C^2 


(a^t-|  +  ct2^' 2 ^  =  ^  G?3  +  G  ^ 


24 


44 


which  can  be  written  in  matrix  form  using  Eqs.  (36)  and  (37)  as 


C11  C12^ 


C21  C22 


f  cf 


F  &I3  +  6  C14 


i  ii.f' 

i  X 


23  ~  G  C24 


The  above  matrix  can  be  easily  inverted  to  give 


and 


G  C-j ^  —^22  +  ^  G23  +  G  G24j1^12 
C11  C22  "  C12  C21 

-F  C23  *  6  C24~*C1 1  *  F  C13  *  6  C14-C21 
C11  C22  "  C12  C21 


These  solutions  can  then  be  substituted  into  Eq.  (35)  to  give  tne 
reflected  field  component  of  the  scattered  field.  The  total  scattered 
field  from  the  flat  plate  is  then  given,  using  Eqs.  (32)  and  (35),  by 


(38) 


Rs(es,=s)  ■  Rd(-S 


•O  +  Rr(*. 


This  solution  has  been  compared  with  measured  results  using  a 
short  dipole  antenna  mounted  above  a  flat  square  plate  as  shown  in 
Fig.  16.  In  Fig.  17  a  dipole  is  mounted  above  the  center  of  a  square 
plate  with  resultant  patterns  shown  for  0°  <_  js  <_  360°  and  es  =  90°. 
In  Fig.  18  the  dipole  is  mounted  near  one  edge  of  the  plate.  In 
Fig.  19  the  geometry  is  the  same  as  the  last  case  except  the  dipole 
is  now  rotated  by  90°.  Note  that  good  agreement  is  obtained  in 
each  case  between  the  calculated  and  measured  results  even  through 
the  back  lobe  region. 

In  order  to  illustrate  the  versatility  of  this  solution,  it  is 
used  to  approximate  the  scattering  effect  of  a  disc.  This  is  done 
by  computing  the  pattern  of  a  monopole  mounted  on  plates  with  in¬ 
creasingly  many  sides.  In  Fig.  20  the  calculated  results  for  plates 
with  4,  6,  8  and  10  sides  are  illustrated  and  compared  with  the 
measured  result  taken  on  the  disk[40].  Note  that  as  the  number  of 
sides  is  increased  the  closer  the  computed  and  measured  results 
agree. 
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Even  though  the  above  results  show  good  agreement  one  must 
realize  the  inherent  approximations  in  this  solution.  It  is  based 
on  edge  diffraction  with  just  singly  diffracted  edge  rays  being 
considered.  Thus,  it  has  been  assumed  that  the  plate  is  large  in 
terms  of  the  wavelength  such  that  double  diffraction  is  normally 
negligible.  However,  neglecting  double  diffraction  may  cause 
some  error  especially  when  the  pattern  is  computed  in  the  plane  of 
the  flat  plate.  Secondly,  a  diffraction  term  from  each  of  the 
corners  should  be  included  but  is  not  available  in  practical  form 
at  present.  Nevertheless,  it  has  little  effect  on  the  overall 
pattern  except  when  the  diffraction  point  approaches  a  corner. 

In  these  two  cases  our  solution  can  be  somewhat  in  error  although 
only  a  small  angular  region  is  involved. 


Fig.  18— Ec  radiation  pattern  for  a  short  dipole  mounted  above  a 

o  n  n 

rectangular  plate  for  e„  =  90°  and  0°  <  <j>  <  360°. 

o  —  S  — 
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Tiis  solution,  of  course,  yields  the  radiation  pattern  for  antennas 
mounted  above  flat  plate  reflectors.  As  stated  earlier  a  large  antenna 
or  array  can  be  handled  simply  by  using  the  superposition  principle  and 
summing  the  scattered  fields  from  delta  function  or  array  elements. 

As  is  shown  in  the  next  chapter,  this  solution  can  be  modified  somewhat 
and  applied  to  the  aircraft  antenna  problem  in  order  to  handle  the 
scattering  by  wings  and/or  horizontal  stabilizers  which  are  approximated 
by  "n"  sided  flat  plates.  Another  possible  application  of  this  study 
is  in  solving  the  near  field  scattering  of  buildings.  This  is  a  problem 
of  recent  interest  in  terms  of  the  overall  airport/aircraft  antenna 
system  performance  in  that  the  scattering  from  buildings  may  adversely 
affect  the  desired  antenna  performance. 

C.  Near  Field  Scattering  by  a  Finite  Cylinder 

Even  though  many  authors^ 41 ,42^  have  previously  studied  the  far 
field  scattering  properties  of  finite  cylinders,  the  near  field  scat¬ 
tering  of  large  finite  cylinders  has  received  little  attention.  It  is 
this  near  field  problem  that  is  examined  here  basically  to  study  the 
engine  effect  in  the  or.-aircraft  problem.  The  geometry  of  this  near 
field  scattering  problem  is  illustrated  in  Fig.  21.  The  solution 
developed  is  basically  an  application  of  ray  optics  techniques. 
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The  nodal  solution  for  the  circular  cylinder  is  used  to  determine 
the  scattering  properties  of  the  cylinder  portion,  wedge  diffraction 
is  applied  lo  the  edges  at  the  end  caps  and  physical  optics  is 
used  to  determine  the  specular  scatter  from  the  end  caps.  In  the 
process  of  developing  the  near  field  scattering  properties,  the  far 
field  problem  is  presented  using  similar  solutions.  Again  this 
far  field  solution  is  obtained  by  applying  several  solutions  which 
are  known  to  be  adequate  in  certain  regions.  Thus,  this  approach 
is  not  an  attempt  to  extend  a  solution  beyond  its  capability  but 
only  uses  previously  developed  solutions  which  are  known  to  be 
valid  in  certain  sectors.  The  most  interesting  point  is  that  there 
is  an  adequate  overlap  between  the  various  regions  in  which  each 
solution  provides  a  good  approximation  thus  giving  complete  coverage. 

Let  us  first  consider  the  physical  optics  solution  applied  to 
determine  the  specular  scatter  from  the  disc-shaped  end  caps.  The 
geometry  used  for  this  analysis  is  illustrated  in  Fig.  22.  In  this 
solution  as  well  as  the  ones  to  follow,  the  parallel  and  perpendicular 
components  of  the  incident  F-field  are  considered  separately.  Con¬ 
sidering  the  perpendicular  component  of  the  incident  field,  one  finds 
that 

-h  a  J*k(z  cos  6.  +  x  sin  e.) 

Ej  =  y  e  1  1 


-r  "j  "  j 

H  =  (x  cos  i  .  -  z  sin  i>.)  j- 


1  jk(z  cosi^.+x  sino^ ) 


The  physical  optics  current  is,  then,  given  by[43] 


jk  x  sin  e.j 

Jx  =  2z  x  H1 |z_0  =  2 y  cos  ^ - 


Fig.  22— Plane  wave  field  incident  on  finite  cylinder  end  cap. 
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ejkp  '[sine-jCos<f  '+sine(cos$cos<t>'+sin4,  si  n4> 1 )  , 


E  =  -2  cos  $  cos  e- 

A  I 


a  r  2v 


-  e 


jkp '[sine^cos^'+sineCcos^  cos<f>'+sin<j>  si r>4> ' l(j 


In  order  to  simplify  this  solution  let 


sin  ^  cos  +  sin- (cos  :  cos  +  sin  :  sin  :') 
=  -  (cos  5  cos  «■  sin  ;  sin  cosU-:') 


where 


J2  2 

sin  9.  +  2  sin  -  sin  e  cos  :  +  sin  -  . 
Substituting  the  above  expression  into  Eq.  (39)  one  finds  that 


Es.  =  -2  cos  9  sin  :  cos  9.  *  ^  eJ*k*'v  cos("'; ' ^dt'dc ' 
x”  1  -0  -0 


Es  =  -4-  cos  9  sin  :  cos  9-  !a  c'  J  (k; 1  % )dc 1 

J.6  TO' 

'0 


or 


(42)  - 


4^a  cos  9  cos  sin  ; 


( ka  v ) 


and  similarly 


(43) 


JL<t- 


4tt3  COS  9.  COS  i 


J-j(kav) 


Let  us  now  consider  the  parallel  component  of  the  incident  field, 
which  is  given  by 


„  jk(z  cos  e.  +  x  sin  e.) 

E?,  -  e  e 


and 


jk(z  cos  e.  +  x  sin  e.) 


"l  =  -y  - 


Again  the  physical  optics  current  is  simply  given  by 


J-  =  2  z  x  H], 


„  2  Jkp1  sin  e-j  cos  <J>‘ 


|  =  x  y~  e 

'  z=0  Lo 


56 


with  the  far  scattered  fields  given  by 


(44)  E.  =■  -2  cos  cos  ;  !  .  '  e 

II  ft  '  _  i  _ 


^a  r2^  jk. ’[  sin  .jCOs,- '+sin  cos( . 


0  '0 

j 

?'  e 


ra  r 2-  jk^,rsin-?.cos:,+siP-cos(i-?,)j 
(45)  E*  =  2  sin  ;  |  1  -  1 

-  0  *  0 


d:'dc' 


These  integrals  can  be  evaluated  using  the  same  technique  as  applied 
earlier  giving 


(46) 
and 

(47) 


rS  _  4-a  cos  •  cos  s 

ii  e  kv 


J-j  (ka-, )  , 


rs  _  4tt3  sin  *  ,  x 

V  "  — -  Jl(Ka,)< 


In  order  to  consider  the  total  scattered  field,  one  must  treat 
an  arbitrary  plane  wave  incident  on  the  end  cap.  In  so  doing  let  us 
assume  that  the  incident  field  is,  in  general,  given  in  ray  form  by 


<48>  ^at  origin  ■»!  i  + 


Then  the  total  scattered  far  field  is  simply  given  using  the  super¬ 
position  principle  and  Eqs.  (42),  (43),  (46),  (47), (48)  by 

(49>  ^end  rW>  '  »<Ri  Efo  *  RX9>  + 

+  *<R1EI«  +  RX>‘ 


Note  that  the  phase  of  the  above  solution  is  referenced  to  the  center 
of  the  end  cap  disc.  If  the  above  solution  was  for  the  rear  end  cap 
with  its  phase  referenced  to  the  center  of  the  cylinder,  then  Eq. 

(49)  must  be  multiplied  by  the  additional  phase  factor 

jk  sl/2  ( cos  e.  +  cos  e) 

Ce  1  ]  . 


57 


It  is  well-known  that  the  physical  optics  solution  is  a  very 
good  approximation  in  the  specular  region  of  the  scattered  field. 

The  center  of  this  region  is  defined  by  setting  Eq.  (41)  equal  to 
zero  (>=0).  This  region  is  bounded  by  the  zeros  of  J](ka. )  as 
discussed  in  Ref.  (44).  However,  these  bounds  may  change  depending 
on  the  size  of  the  disc  which  in  turn  varies  the  region  of  overlap 
between  the  physical  optics  solution  for  the  end  caps  and  the  wedge 
diffraction  solution  for  the  circular  edges  of  the  finite  cylinder. 
In  any  event,  the  validity  of  this  approach  is  illustrated  by  the 
results  presentee  later. 

Now  the  question  arises  as  to  how  this  solution  can  be  applied 
to  the  near  field  scattering  problem.  Note  that  this  near  field 
problem  will  ultimately  be  adapted  to  simulate  the  scattering  from 
aircraft  engines.  Even  though  this  flat  structure  does  not  simulate 
the  ends  of  engines,  it  is  a  convenient  means  of  closing  off  the 
surface.  Further,  as  shown  later  the  engine, has  very  little  effect 
on  the  over-all  pattern.  In  general,  aircraft  type  engines  are 
inherently  long  in  comparison  to  their  diameter;  thus,  the  source 
can  be  in  the  near  field  of  the  engine  yet  in  the  far  field  of  the 
relatively  small  discs  at  the  ends  of  the  engine.  As  a  result,  the 
far  field  solution  of  Eq.  (49)  can  be  adapted  to  this  problem  simply 
by  assuming  that  field  intensity  of  the  incident  plane  wave  is  that 
radiated  by  the  source  and  incident  at  the  center  of  the  end  cap. 
This  simply  implies  that  the  field  incident  at  any  point  across  the 
end  cap  is  very  nearly  the  same  as  the  field  at  the  center  of  the 
end  cap  (or  nearly  a  plane  wave).  The  geometry  for  this  situation 
is  illustrated  in  Fig.  23  with  the  source  located  behind  the 
cylinder. 


As  was  done  for  the  flat  plate  problem  it  is  known  that  the  far 
field  pattern  of  a  general  source  can  be  wntten  as 


^sojrce^a’ ^'Ja  ^a*'^  +  ^a^ea,<J>a^ 


jke 


-jkr 


4nr 


where 


jke 


-jkr 


4nr 


is  factored  out  as  done  earlier.  This  field  is  then 


decomposed  into  parallel  and  perpendicular  components  such  that 


(50) 


fsource<V*a>  =  &<Fx+  M 


which  is  the  field  incident  on  the  rear  end  cap.  The  vr.lues  of  ea  and 
<t>a  are  known  once  the  center  of  the  end  cap  is  located  in  terms  of 
the  source  coordinate  system.  The  total  near  field  scatter  from  the 
rear  end  cap  is  then  given  using  Eqs.  (42),  (43),  (46),  (47),  (50)  by 
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Fig.  23— Near  field  specular  scatter  by  rear  cylinder  end  cap. 


(51) 


xe 

jkR 


Jrac-  « 


3jk  %! 2  cos  ei 


for  kay  £  (approximately)  3.8317.  Note  that  in  this  case  the  phase 
is  referenced  to  the  center  of  the  cylinder.  A  similar  expression 
can  be  written  when  the  source  is  located  in  front  of  the  cylinder 
where  it  is  possible  that  the  specular  scatter  from  the  front  disc 
can  dominate. 


It  is  usually  extremely  difficult  to  analyze  the  more  complex 
end  shapes  of  actual  engines.  Thus,  one  must  have  a  real  need  for 
such  results  to  make  a  thorough  investigation  of  that  topic.  As  will 
be  shown  in  the  following  chapter  the  engine  has  little  effect  on 
the  principal  plane  patterns.  Consequently,  the  flat  plate  end  cap 
simply  serves  to  close  off  the  structure. 


Once  removed  from  the  specular  region  of  the  end  caps,  the 
diffractions  from  the  circular  edges  become  dominant.  In  order  to 
include  these  effects,  the  problem  requires  more  sophistication  than 
the  striaght- forward  wedge  diffraction  problem  of  Section  II-C. 

The  curvature  of  the  diffracting  edges  must  be  included  in  this 
solution  by  applying  an  equivalent  current  technique  such  as  that 
presented  in  Ref.  j_45I.  Using  this  approach  an  equivalent  current 
is  found  on  the  diffracting  edge  using  the  appropriate  two- 
dimensional  diffraction  solution.  The  resulting  current  is  then 
used  in  the  radiation  integrals  to  find  the  far  scattered  field. 


Let  us  first  consider  the  perpendicular  component  of  the 
incident  plane  wave  on  the  back  edge  of  the  cylinder  as  illustrated 
in  Fig.  22.  This  component  is  given  by 


E1  = 


jk(z  cos  e.  + 

y  e 


x  sin 


which  results  in  giving 


z=0 


cos  *'  e 


jkx  sin  e- 


and 


(53)  h] 


z=0 


cos  e.  sin  jkx  sin 

7  -  -  -  0 


The  equivalent  electric  and  magnetic  currents  are  given,  respectively, 
by  Ref.  [20]  as 


(54) 


(55) 


Ie  ,  =  E1 

x*  Z  k  sin  B  x* 


7=0  ’ 


and 


I 


m 

V 


2jZGh4')  , 


- ~T~ 

k  sin  b 


H' 


|z=0  • 

which  are  related  to  the  ray  form  of  the  diffraction  function.  The 
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various  terms  of  the  above  expression  are  defined  later  in  the  dis¬ 
cussion.  In  any  event,  these  currents  are  qiven  using  Eqs.  (52)-(55) 
by 


=  2j  Ge(;')cos 
2ol  si"2= 


jka  sin  •  .  cos 
e  1  ,  and 


i".. 

j.  ; 


2jGh(i 1 )  cos 


sin 


k  sin  3 


jka  sin  cos  c* 
e 


Substituting  these  currents  into  the  radiation  integrals  one  finds 
that 


jka[sine.cosj  '+sinecosi(i-i')] 


r  2tt 

Ee.  =  -a  cos  0  1  Je  ,  sin(^')e 

16  Jo  ■L? 

{27  jka[sine.cos; '+sin6cos(i-i' )j 

E®  =  -a  J®,  cos(;-;')e  1  do’ 

■LQ  1  n  V 


d;V 


(56) 

EmA  =  a  cos  6 


27 


Jm  .  sin($-*')e 


jka[sine.cos*  '+sinecosv$-$')] 


J0  ^ 


do' 


^9  *  f0 


jka[sine.cos<*>'+sinecos(<+>-<j' )] 


do' 


where 

(57) 


oe  (,.j  =  and 

k  sin^B 

m  2j  Gh(<*>')  cos  e.  sin  o' 

J>')= . - 


k  sin  8 


Note  that  the  phase  reference  is  at  the  center  of  the  end  cap  and 

the  factor  e~'^*<rJ  has  been  removed  from  the  above  expressions. 

Now  using  the  method  of  stationary  phase  one  can  evaluate  the  previous 
integrals  asymptotically  usingC46] 


61 


!  F(:,)eJKa’'’  jd;' 
J0 


F(;;) 


2- 


ka|  :■“(;')  | 


j!  ka  ^ ( : ^ )  +  ^  sin  ?"(  ^)  J 

.  e 


where  is  the  stationary  phase  point  or  points  defined  by  :'(vs)=0- 
All  of  the  above  integrals  { Eq .  (56))  are  of  this  form  with 


*(;')  =  sin  e.  cos  j'  +  sin  e  cos(j-$') 
giving  the  relation  that 


tan 


sin  v  sin  ; 


vSin  +  sin  A  cos 


which  has  two  solutions  for  0  ^  <  360  ( *51 ,  $52)*  Making  the 

appropriate  substitutions  in  the  above  integral  expressions  one  finds 
that 


(58)  Ee  (0,0)  -v  -a  cos  9  \  JeA .  Uo)sinU-$£)  X(^) 

<1*7  a  x  -1-  V 

:SP'S2 


E*  (»,♦) 


-a  l  ^.(♦s)c°s (♦-*$) 

*S1’*S2 


and 


Em(0,$)  ^  a  cos  9 

-L<P 


l  O^.UpsinU-^)  X( 4>g) 
*S1**S2 


E^0 (0 ,4>)  ^  -a  1  J^i  (♦£)  cosU-ifrg)  ^Uj) 


PS'I5?S2 


J-<P 
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XU 


*’ 7 
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jCka*(*g)  +  J  sign  *"US)] 


ka|*"(*q)| 


62 


and 


:(vo)  =  -$" ( tc)  =  sin  e.  cos  +  sin  0  cos(;-t^-). 


The  total  perpendicular  component  of  the  scattered  far  field  for 
the  circular  diffracting  edge  is  then  the  sum  of  the  electric  and 
magnetic  components  and  is  given  using  Eqs.  (57)  and  (58)  by 


(59)  E^,;)  =  2ja  > 


;S1’*S'  k  sin  3S 


x  ( *,  • ) 

-^2 — [-Ge(^)cose  cos^in(:-:£)  + 


L 

G  (;£)  cos  sin  cos(j-c^)] 


d  x(^c) 

(60)  E^(6,*)  =-2ja  l  —^2 


»S1’«S2  k  s1n  6S 


[G  (a^)cos  ^  cos(*-^)  + 


I* 

Gn(^)cos  e  cos  91-  sin  $1  sin($-<^)]  . 


The  diffraction  terms  are  yet  to  be  defined  in  terms  of  the 
stationary  phase  points  ( <+>si » '+>$2) •  The  geometry  pertinent  to  this 
discussion  is  illustrated  in  Fig.  24.  The  ray  forms  of  the  dif¬ 
fracted  field  are  given  by 

e 

Gh(^)  =  R(*-<|-0,3/2)  +  RU+*q,  3/2) 


where  it  is  shown  in  Ref.  [20]  that 
|  sin  — ■ 

R<^2>  =V - 25 

cos  2 —  cos  y*- 

The  values  of  the  diffraction  angles  are  based  on  the  two-dimensional 
wedge  problem  and  are  given  by 


*0  =  tan 


-1 


r-I  •  y 


A 

<-I  •  X. 


el  =  tan"1 


COS  0. 


.-cos  sin  e7.  / 
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where  0  1  180°.  Fina11y>  if  tne  condition  exists  where  .  or 
jq  >  n-,  then  diffractions  do  not  occur  from  '"hat  stationary  phase 
point  in  that  the  cylinder  shadows  tnat  terr- 


The  parallel  component  of  the  incident  field  is  given  by 


.  „  jk(z  cos  *•  +  x  sin  e^) 

E?,  =  (x  cos  -  z  sin  e.)e  ! 


which  gives 
-i  i 

‘'z=0 


jkx  sin  *<- 

E ‘  i  =  -cos  6-  sin  e  ,  and 

ilvU_r  1 


H' 


•'6'z=0 


COS  C  1 


jkx  sin 
e 


Then  using  the  same  approach  as  before,  the  parallel  component  of 
the  scattered  field  for  the  diffracting  edge  is  given  by 


(61)  E.d  >,,)  =  2ja  l 


X(:S>  „e 


ne 


"sii0S2  k  sin  eS 


[G  (r.^)cos  e,  sin  sin(6-;j) 


and 

(62) 


+  Gh(<^)cos  cos(>-tj)] 


X(^) 


Ej  (o,*)  =  2ja  l - [GeU£)cos  9i  sin  cos($-$£) 


^S1 ’^S2  k  S1P  6S 


-GnU^)cos  e  cos  ^  sinU-^)] 


with  the  various  terms  being  defined  previously. 

The  total  scattered  far  field  from  the  rear  diffracting  edge 
is  given  using  Eqs.  (48),  (59),  (60),  (61),  (62)  by 


65 


(63) 


rfa r  field  / 
edge  diff . ' 


;  v-}i  *  v> + ; 

<*V  .  R*  Cd  ,  'e'^  ^ 

j.  j.  ;  it  11  - 


Note  that  in  this  case  the  phase  is  referenced  to  the  center  of  the 
cylinder.  In  addition,  a  similar  expression  can  be  generated  for 
the  scattered  field  which  results  from  diffractions  oft  the  front 
edge  of  the  finite  cylinder.  These  two  diffracted  field  solutions 
are  then  summed  using  the  superposition  principle  and  applied 
between  the  specular  regions  of  the  end  caps  and  cylinder  section. 

The  cylinder  specular  term  is  presented  later  in  the  discussion. 

It  was  assumed  by  using  the  static  .ary  phase  method  that  only 
two  points  around  each  of  the  cylindrical  rims  could  contribute  to 
the  scattered  field.  However,  there  are  certain  regions  (caustics) 
where  diffractions  from  the  complete  rim  contribute  to  the  scattered 
field.  In  these  regions  one  must  carry  out  the  complete  integrations 
around  the  edges  rather  than  using  the  stationary  phase  method. 

The  expressions  have  been  developed  but  are  not  included  in  our 
numerical  computations  in  that  this  solution  is  applied  only  to 
simulate  the  engine  near  field  scattering  effect.  Thus,  a  complete 
investigation  of  the  finite  cylinder  is  beyond  the  scope  of  our 
final  objective.  For  this  reason,  our  numerical  solutions  can  tend 
to  diverge  near  the  caustic  regions  as  is  noted  when  the  results  are 
presented  later. 

As  was  done  for  the  end  cap  specular  region,  it  will  be  assumed 
that  the  edges  of  the  finite  cylinder  are  in  the  far  field  of  the 
source  for  the  near  field  scattering  problem.  With  this  being 
the  case  and  using  the  geometry  illustrated  in  Fig.  22,  the  dif¬ 
fracted  field  from  the  rear  edge  is  given  using  Eqs.  (50),  (59),  (60;, 
(61),  (62)  by 


(64)  I3E  =(<<F Me  *  WO  + 

-jkFt 


C(F1+GJE“t  +  (F.+GX^  ^ 


jke 


Jk  (ft/2)  cos  e 


A  similar  expression  can  be  found  for  the  diffracted  field  from  the 
front  edge  of  the  cylinder.  The  total  diffracted  field  is  then  the 
superposition  of  the  diffracted  fields  from  the  two  edges  of  the 

r\/l  t nHov'. 

V  "  * - -  ■  • 
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Since  in  this  case  the  range  of  the  source  is  not  an  infinite 
distance  from  the  edge,  one  must  apply  the  form  of  the  diffraction 
function  in  which 


R(:  ,3/2) 


W 


.  2 
Sin  r 


S’ 


:  ,3/2 )e 


j(h.Rrsin 


*  nI  2-kR 


sin  r., 


giving 


j(kRt.  sin2^ 
e 


+-/4)  _ 

J 2-kRr 


sin 


S* 


One  then  simply  subst  tutes  this  result  for  the  terms  EL  _ ,  E^.,  Eir, 
and  E$-  for  the  near  ^ield  scattering  solution  of  Eq.  (64).  Otherwise 
the  terms  are  the  same  as  found  previously.  Finally,  in  the  near  field 
problem  the  diffracted  fields  from  the  two  einite  cylinder  edges  are 
included  everywhere  except  in  the  specular  ^egion  of  the  end  caps. 

This  is  somewhat  different  from  that  of  the  far  field  solution  v  .ere 
the  diffracted  Melds  are  only  present  between  the  specular  regions 
of  the  end  cap  and  cylindrical  sections. 


The  final  region  to  be  considered  is  th''  specular  region  of  the 
long  cylindrical  section.  The  solution  for  this  region  is  found 
similar  to  the  physical  optics  solution  used  for  the  end  caps; 
however  in  this  case,  the  current  used  in  the  radiation  integral  is 
the  exact  current  found  on  an  infinitely  long  cylinder.  This  current 
is,  then,  integrated  over  the  finite  length  of  the  cylinder  as 
presented  in  Ref.  [42].  This  type  of  solution  generally  has  good 
accuracy  throughout  the  specular  region  of  the  cylindrical  section  much 
as  the  physical  optics  solution. 

With  a  perpendicular  component  E-field  incident  on  the  cylinder, 
it  is  shown  in  Appendix  I  that 
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where  ■  =  k  sin  ..  The  equivalent  surface  current  density  is  qiven 
by  '  1 

0  =  r  X  (H*  z  +  : )  |  =  (.=a)  +  z  H^(.=a) 

jl  J.Z  j.4*  1  _«  j-z  -- 

The  component  of  the  far  scattered  field  is  given  by 
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where  Xn(a)  =  Jn(na)  -  -T^y 
evaluated  using  uK  ' 
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and 


ejka  Si"  =  2-ej"(:+-/2)  j  (ka  sjn 


Substituting  those  results  into  Eq.  (65)  one  finds  that 
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Similarly  the  .--component  of  the  scattered  field  is  given  by 
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which  can  be  evaluated  as  above  to  give 
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For  a  parallel  component  E-field  incident  it  is  shown  in 
Appendix  I  that 


H  =  0 

llZ 

f  1  j  k2  COS  6.  00  J  (Ka) 

H...  --  4-  e  l  jn  CJ^(Bp)  -  4or—  h!2)  (£p):e“jn  . 


'<■*  l „ 


n=-co 


Hp^(8a)  n 


69 


In  this  case,  the  equivalent  surface  current  density  is  qiven  by 
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Substituting  this  current  into  the  radiation  integrals  one  finds 
that 
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This  integral  can  De  evaluated  as  done  previously  giving 
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The  total  scattered  field  in  this  region  is  then  given  using  the 
superposition  principle  and  Eqs.  (66),  (67),  (69),  (70)  by 
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where  Rx  and  R„  are  defined  in  Eq.  (48).  Note  that  in  this  case 
the  phase  is  already  referenced  to  the  center  of  the  finite  cylinder. 


As  stated  earlier  this  solution  gives  a  good  approximate  solution 
for  computing  the  specular  scattered  fields  of  the  long  cylindrical 
section.  As  was  done  earlier  this  region  is  normally  defined 
between  the  first  nulls  about  the  main  beam  of  the  specular  scatter 
directior.  These  boundaries  are  each  defined  by 

t£-  (cos  e.  +  cos  e)  =  nn 
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where  n=0  gives  the  SDecular  main  beam  direction  and  n  =  *1  defines  the 
first  null  directions.  This  is  the  only  term  necessary  through  the 
cylindrical  specular  scatter  region.  However,  as  is  shown  later  the 
edge  diffracted  fields  must  be  included  through  this  region  for  the 
near  scattered  ^ield  solution. 


This  solution  completes  our  study  of  the  scattered  far  field 
from  a  finite  cylinder  which  is  given  using  Eqs.  (49),  (63),  (71)  by 
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In  order  to  compare  this  solution  with  previously  measured  scattered 
far  field  data  it  was  found  that  backscatter  data  was  most  readily 
available.  In  Fig.  25  the  results  of  Eq.  (72)  are  compared 
with  measured  bacxscatter  data  taken  from  Ref.  [413.  It  was 
extremely  difficult  to  extract  the  measured  data  from  their  graphs 
in  that  the  finite  cylinder  measured  was  quite  large  which  resulted 
in  much  ripple.  In  any  event,  the  agreement  between  the  results  for 
the  points  chosen  does  indicate  the  validity  of  our  solution.  Note 
that  both  incident  polarizations  were  considered  which  covers  the 
complete  problem  at  least  for  the  backscatter  case. 


The  far  field  solutions  could  be  applied  to  the  near-field 
problem  in  the  two  other  regions  since  the  cylinder  was  assumed  to 
be  long  and  thin.  However,  that  approximation  can  not  be  made  for 
this  cylindrical  specular  region  unless  the  source  is  moved  a  large 
distance  (2;:2/x)  away  from  the  cylinder.  This  implies  that  the 
spherical  wave  front  radiated  by  the  source  must  be  considered  in 
the  approximating  currents.  This  is  a  very  difficult  problem  to 
solve  except  for  cylinders  short  in  terms  of  the  wavelength.  Thus, 
a  different  approach  must  be  applied  to  solve  for  the  scattered  near 
field  solution  in  this  specular  region. 

This  obstacle  is  overcome  by  applying  a  technique  similar  to 
that  employed  for  the  finite  flat  plate  problem,  in  which  the  re¬ 
flections  from  the  plate  are  added  only  when  the  point  of  reflection 
falls  within  the  limits  of  the  finite  plate.  In  this  case,  the 
modal  solution  can  be  used  to  obtain  a  similar  scattered  field  term 
from  an  infinite  cylinder  which  is  only  included  when  the  c^gin  of 
the  scattered  term  falls  within  the  limits  of  the  finite  cylinder. 


Based  on  a  ray  optics  description,  it  is  known  that  for  the  reflecteo 
terms  the  angle  of  incidence  (••■,-)  is  equal  to  the  angle  of  reflection 
(f).  Secondly,  it  is  known  that  the  creeping  wave  terms  attach  to 
the  cylindrical  surface  tangentially  and  propagate  around  the  surface 
in  a  helical  path  with  constant  pitch,  such  that  •>  remains  constant. 
For  these  reasons,  one  can  affirm  that  to  a  good  approximation  the 
bounds  of  this  region  can  be  specified  in  terms  of  as  shown  in 
Fig.  26,  for  :•  1  -  -  l  ~f. 


LOCATION 

Fig.  26— Cylindrical  specular  region  for  near  zone  cylinder. 


In  order  to  apply  the  results  of  Appendix  I  one  must  use  a 
reciprocity  argument  in  that  it  is  assumed  there  that  a  plane  wave 
is  incident.  In  this  case,  the  source  is  in  the  near  field  and 
the  scattered  far  field  is  desired  which  can  be  solved  using 
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where  8=k  sin  «■;.  Note  that  the  above  equations  give  the  9  and  $  com¬ 
ponents  of  the  field  at  the  source  position  which  is  defined  by  (p,$,z). 

The  received  signal  by  the  antenna  can  be  found  using  the  far 
field  pattern  as  defined  by 
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Thus,  the  pattern  function  of  the  source  can  be  expressed  in  terms 
of  the  e  and  $  coordinates  by 
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where  ea  and  $a  are  defined  in  terms  of  the  incident  field  direction 
in  the  source  coordinate  system.  The  radiated  field  by  tne  source 
using  the  reciprocity  theorem  is  then  given  by 
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for  -r  ■=  -j  l  f  as  defined  in  Fig.  26  and  where  and  now 
define  the  scattered  field  direction.  In  this  case  the  phase  is 
already  referred  to  the  center  of  the  finite  cylinder.  In  addition, 
the  diffractions  from  the  edges  of  the  cylinder  must  be  included 
throughout  this  region.  Note,  also,  that  the  directly  radiated 
field  by  the  source  must  be  included  in  every  region  of  the  near 
field  problem. 

To  illustrate  the  validity  of  this  near  field  solution,  it 
is  compared  with  measured  results  for  two  cases.  A  vertical  and 
horizontal  dipole  are,  respectively,  mounted  in  the  near  zone  of 
a  finite  cylinder  as  shown  in  Fig.  27a  and  b.  One  should  note 
that  the  agreement  is  quite  good  except  near  the  caustic  regions  where 
our  solution  tends  to  diverge.  The  reason  is  that  the  diffractions 
from  the  circular  edge  are  not  coming  from  just  two  points  but  from 
the  complete  circular  edge,  which  can  be  included  in  our  solution  by 
actually  integrating  Eq.  (56).  However,  it  is  shown  in  Section  IV-D 
that  when  this  solution  is  used  in  solving  for  the  azimuth  principal 
plane  pattern,  the  engine  (finite  cylinder)  has  very  little  effect. 
Thus,  an  extensive  study  of  this  topic  is  not  necessary  at  present. 

In  any  event,  the  basic  tools  are  presented  in  the  previous  discussion. 


Fig.  27a— Short  dipole  illuminating  a  finite  cylinder. 
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Fig.  27b— Radiation  pattern  of  a  horizontal  dipole  (Ej 
with  -180°  <  o  <  180°. 
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V 


CHAPTER  IV 

ON-AIRCRAFT  ANTENNA  PRINCIPAL  PLANE  PATTERN  STUDY 


A.  Introduction 


This  chapter  illustrates  the  application  of  solutions  previously 
discussed  to  determine  radiation  patterns  in  three  principal  Dianes 
for  fuselage  mounted  aircraft  antennas.  Most  of  thest  results  have 
been  presented  elsewhere[47,48l;  nowever,  tr.ey  are  pertinent  to  the 
present  discussion.  Furthermore,  several  improvements  have  been 
made  in  the  present  work.  The  complete  detailed  equations  for  each 
solution  will  not  be  presented  except  in  the  cases  where  improvements 
have  been  made. 

The  basic  aircraft  to  be  analyzed  in  this  study  is  shown  in 
Fig.  28.  It  ’s  composed  of  flat  plates,  cylinders,  cones,  and  spheres. 
It  is  assumed  that  the  source  is  mounted  on  the  fuselage  and  re¬ 
stricted  to  the  regions  near  the  top  or  bottom  of  the  aircraft.  As 
is  shown  later  in  this  discussion  some  of  these  restrictions  and 
models  are  revised  in  order  to  obtain  better  approximate  solutions 
for  the  desired  pattern.  In  this  way  models  are  considered  that 
approximate  a  wide  variety  of  aircraft  structures  with  the  solutions 
derived  in  such  a  form  that  arbitrary  antennas  can  be  considered 
simply  by  integrating  the  equivalent  aperture  currents. 

The  lower  frequency  limit  of  these  solutions  is  dictated  by  the 
ray  optics  format  which  -equires  that  the  various  scattering  bodies 
be  no  closer  than  approximately  a  wavelength  with  the  overall 
aircraft  being  large  in  terms  of  the  wavelength.  The  upper  frequency 
limit  is  dictated  by  the  model  representation  of  the  actual  aircraft 
considered. 

B.  Roll  Plane  Analysis 

This  section  is  a  synopsis  of  the  material  presented  in  Ref.  [47 j 
along  with  recently  developed  improvements.  Tne  two-dimensional 
problem  is  considered  initially  in  order  to  develop  the  necessary 
analytical  tools  to  attack  the  much  more  difficult  three-dimensional 
roll  plane  problem.  The  geometry  of  the  two-dimensional  problem  is 
illustrated  in  Fig.  29  with  and  without  the  engines  included.  The 
fuselage  and  engines  are  assumed  circular  in  cross-section  and 
mounted  symmetrically  aout  the  finite  wing.  Since  wedge  diffraction 
is  applied  to  handle  the  finite  length  wing,  the  radiated  field  must 
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STATION 
1-2  SPHERE 

-3  CIRCULAR  CYLINDER 
-4  CONE 

5  CIRCULAR  CYLINDER 

6  PLANAR  CONDUCTOR 

7  PLANAR  CONDUCTOR 


Fiq.  28— Simplified  aircraft  model 


be  described  in  terns  of  rays.  However,  one  of  the  nicer  features  of 
this  approach  is  that  other  solutions  such  as  nodal  solutions  can  be 
cast  into  a  ray  fcm  and  then  applied  to  a  diffraction  problem.  Con¬ 
sequently,  it  was  found  that  the  nodal  solutions  (Eos.  (3)  and  (4))  for 
an  arbitrary  antenna  counted  or,  an  infinite  circular  cylinder  was  ideal 
for  treating  the  antenna  counted  cn  the  fuselage.  In  fact,  this 
solution  has  been  applied  in  past  years  as  the  sole  solution  for  liigr 
frequency  on-aircraft  antenna  analyses. 

Using  image  theory  the  wing  reflections  are  analyzed  using  the 
geometry  illustrated  in  Fig.  30.  liote  that  due  to  the  finite  length 
wing  the  reflected  field  contributes  only  in  a  given  region.  The 
wedge  diffraction  solution  is  applied  to  include  the  effect  of  the 
wing  tip.  The  field  scattered  by  the  engines  is  included  using  the 
modal  solutions  {Appendix  I)  for  a  near  field  source  illuminating  an 
infinitely  long  circular  cylinder.  In  addition,  a  portion  of  the 
energy  scattered  by  the  engine  is  incident  upon  the  wing  tip  which 
is  included  as  a  higher-order  tern.  The  phase  of  these  various  terms 
is  then  referenced  to  the  center  of  the  fuselage,  and  the  terms  summed 
to  give  the  total  radiated  field.  There  were  other  higher-order  terms 
considered  in  this  study  but  found  to  be  negligible  so  they  are  not 
included  here  or  in  the  following  numerical  results.  There  are 
basically  three  infinitesimal  sources  considered  in  this  analysis 
{monopole  and  axial  and  circumferential  slots).  These  solutions 
allow  for  an  arbitrary  antenna  to  be  considered  using  a  numerical 
aoerture  integration  as  presented  in  Section  II-B. 

The  radiation  patterns  for  three  fini.e  sources  are  shown  in 
Fig.  31  on  a  model  without  engines.  The  monopole  is  approxi¬ 
mately  >./4  and  the  slots  are  simply  open-ended  X-band  waveguides. 

These  solutions  compare  very  favorably  with  measured  results  taken 
on  the  two-dimensional  aircraft  model  of  Fig.  29.  It  is  observed 
that  the  direct  term  from  the  antenna  has  the  dominant  effect  on 
the  pattern  in  the  lit  region.  The  reflected  term  adds  the  slowly 
varying  ripple  to  the  lit  portion  of  the  pattern.  The  scattering 
from  the  wing  tips  causes  the  rapidly  varying  backlobes  in  the 
pattern.  The  backlobes  are  not  shown  in  these  figures  in  that 
they  were  well  below  the  noise  level  of  our  equipment  and  thus,  the 
measured  result  was  not  accurate  through  that  region.  The  radiation 
patterns  for  the  same  configurations  with  the  engines  added  are 
illustrated  in  Fig.  32.  Again  very  good  agreement  is  obtained 
between  our  calculated  and  measured  results.  These  results  indicate 
that  the  scattering  from  the  engines  tend  to  smooth  out  the  pattern  by 
filling  the  nulls  and  lowering  the  peaks  especially  near  =  90°  and 
270°.  ' 
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MEASURED 


radiation  pattern  for  a  x/4  monopole  on  a  fuselage 
i  dominant  engine  and  finite  length  wing. 


Even  though  the  two-dimensional  solution  should  be  more  accurate 
in  some  cases  than  the  previously  used  modal  solutions  which  only 
models  the  fuselage,  one  must  model  the  wing  width  as  well  as  the 
length  in  order  to  make  the  solution  really  practical.  In  order  to 
accomplish  this  feat,  the  near  field  flat  plate  scattering  solution  of 
Section  III-B  is  adapted  to  this  new  model  such  as  illustrated  in 
Fig.  28.  Note  that  each  wing  can  be  located  arbitrarily  with  any 
number  of  edges  provided  only  that  the  wings  are  flat  and  horizontal. 

The  following  discussion  is  presented  with  more  detail  in  that 
selected  improvements  have  been  made  over  our  previous  solution  as 
presented  in  Ref.  G47D. 

Let  us  first  find  the  effective  source  location  for  the  reflected 
field.  Recall  that  in  our  flat  plate  result  the  source  was  imaged 
and  the  reflected  field  added  to  the  total  solution  provided  the 
image  ray  passed  through  the  finite  flat  plate  (wing)  limits.  So 
one  must  initially  determine  the  effective  source  position  which  in 
turn  results  in  giving  the  desired  reflected  field.  With  the  source 
mounted  on  an  infinitely  long  circular  cylinder  one  can  easily  show 
that  the  surface  rays  from  the  source  propagate  around  the  cylinder 
on  helical  paths  (geodesics),  which  in  turn  diffract  energy  tangentially. 
Now  let  us  assume  that  the  source  does  not  illuminate  the  right  wing 
directly  (as  illustrated  in  Fig.  33)  and  proceed  to  determine  the 
unique  helical  path  that  diffracts  energy  from  a  known  tangent  point 
which  is  then  reflected  off  the  wing  in  the  desired  radiation  (or 
scatter)  direction.  This  helical  path  can  be  specified  in  general 
by 


x  =  af  cos  $ 

y  =  af  sin  $ 

z  =  b(<f>-d>  )  +  z 
vso'  sc 


where  (af,  eso,  zSq)  defines  the  actual  source  location  on  the  fuselage. 
Let  the  desired  radiation  direction  be  given  by  (es,<j>s).  Then  using 
Fig.  33a  one  finds  that  the  tangent  direction  at  tne  effective  source 
location  is  defined  by 
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Fig.  33b— Reflection  problem  in  y-z  plane 
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which  result  in  $e  =  90°  -  $  .  Using  Fig.  33b  the  tangent  direction 
is  given  by  s 
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which  gives  b  =  a*  cos  4>e  tan  (90°-es).  The  above  equations  uniquely 
specify  the  effective  source  position  for  reflet;tions  from  the  right 
wing  in  terms  of  the  given  radiation  direct  .‘on  by 

xe  =  af  cos (90°  -  <hs) 

(76)  y0  =  af  sin(90°  -  <f>s) 

ze  =  af  cos (90°  -  <ts)tan(90°-es)('.f/2-(j.s-  gQ)  +  zg0 


which  in  turn  can  be  used  in  the  flat  plate  problem  as  the  effective 
source  location.  Note  that  as  the  desired  radiation  direction  is 
varied  the  effective  source  location  changes.  In  addition,  if  the 
source  directly  illuminates  the  wing  for  a  given  reflection  term  then 
the  effective  source  location  is  simply  the  actual  source  location. 

A  result  similar  to  Eq.  (76)  can  be  found  for  the  reflections  from 
the  left  wing.  Finally,  the  actual  source  field  value  used  to 
compute  the  reflected  term  is  determined  from  the  modal  solutions  of 
Appendix  I. 

Using  a  similar  technique  the  effective  source  locations  for  the 
diffracted  'ield  components  must  be  found.  Recall  that  our  flat  plate 
solution  used  a  search  technique  to  find  the  diffraction  point  by 
computing  the  diffraction  angles  at  selected  test  points  along  a  given 
edge.  Once  a  test  point  (x^,  yd,  Zrf)  is  specified  along  the  edge  one 
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can  find  the  effective  source  location  [y.Pt  ye,  ze)  using  the  geometry 
illustrated  in  Fig.  34.  Again  it  is  assurer  that  the  source  does  not 
directly  illuminate  the  test  point.  Using  F^g.  34a  one  finds  that  the 
tangent  direction  at  the  effective  source  is  given  by 


dx 


(xp>y*»zJ 


-af  sin  *  =  -y 


d* 


(xe*ye,^e) 


=  cos  *e  =  xe 


or 


dx 

2y|/ 


_  <xd  -  se|  _  ^ 

ixe,ye,ze)  X| 

2  .  2  .  .2 


whicn  gives  ^  =  x“  +  ye  and  af  =  XgX^  +  yeyd*  fhe 
(b)  of  the  helical  path  on  the  cylindrical  surface  i 
Fig.  34b  by 


value  of 
s  defined  using 


dz 

<h 


(xe>ye>zp) 


=  b 


or 
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dz 
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"  af  cos  =  xe 
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e 

that  for  the  wing  on  the  right  side  the  effective  sour  ms  location  for 


which  gives  b  =  xQ  (  ]  .  Using  the  aho’-e  relations  one  finds 

side  the 

the  test  diffraction  point  is  given  by 
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\TEST  DIFFRACTION 
POINT 


Fig,  34a— Diffraction  problem  in  x-y  plane. 


Fig,  34b— Diffraction  problem  in  y-z  plane. 


x*  =  tfXd  +  yd^d  +  Xd ' a* 

(77>  j’e  ”  Jx2  J  ~7L j  [afyd  '  xd'i  5d  +  yd  “  af 

7  Vd  (V;so)  *  "so^d-ye1 
e  xel  W  +  lW 

-1  /y.\ 

where  $e  =  tan  )  *  G^ven  effective  source  location  for  the 

chosen  test  point,  the  search  technique  of  Section  III-B  is  applied 
to  find  the  actual  diffraction  point  along  a  given  edge,  h’ote  that 
once  the  actual  diffraction  point  is  determined,  che  effective  source 
of  the  diffracted  field  is  specified  by  Eq.  (77),  and  the  source 
field  value  is,  again,  computed  using  the  modal  solutions. 

The  total  field  is  found  by  summing  the  directly  radiated  field 
with  the  scattered  fields  from  the  wings  using  the  superposition 
principle.  Several  different  configurations  have  been  tested  using 
this  solution  and  compared  with  measured  results.  The  roll  plane 
radiation  patterns  for  the  three  infinitesimal  sources  are  shown 
in  Fig.  35,  when  the  antennas  are  mounted  directly  above  the 
wings.  The  patterns  are  shown  in  Fig.  36  when  the  antennas 
are  mounted  on  the  fuselage  over  the  back  limit  of  the  wing.  In 
Fig.  37  the  patterns  are  shown  when  the  antennas  are  mounted 
over  the  wings  but  rotated  45°  from  the  straight  up  direction.  Only 
the  monopole  is  considered  for  measured  comparisons  since  it  was 
the  only  source  available  at  the  time.  In  each  of  these  cases, 
one  can  note  the  good  agreement  obtained  between  the  calculated  and 
measured  results.  The  slight  deviations  being  attributed  to  the 
assumptions  made  in  solving  the  flat  plate  problem  of  Section  III-B. 
For  example,  the  jump  in  the  calculated  result  of  Fig.  36  near  =  90° 
and  270°  can  be  attributed  to  the  lack  of  doubly  diffracted  terms  in 
the  solution. 

It  was  originally  assumed  in  the  aircraft  model  of  Fig.  28  that 
the  wings  ware  attached  such  that  they  were  contained  in  the  horizontal 
plane  which  also  contains  the  axis  of  the  fuselage.  This  meant  that 
no  diffractions  occurred  from  the  junction  edge  in  that  image  theory 
could  be  applied  to  handle  the  cylinder/plate  junction  problem.  This 
is  not  the  case  when  the  wing  is  moved  away  from  this  location. 
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However,  in  order  to  make  th’s  solution  even  mo'e  practical  one  must 
be  able  to  move  the  wing  from  its  former  centra''  position.  As  a 
result,  a  technique  has  been  developed  to  handle  the  case  with  the 
wing  mounted  above  the  central  position  or  closer  to  the  source 
location.  Again  the  flat  plate  analysis  of  Section  III-B  is  applied 
but  in  this  case  the  edge  formed  by  the  junction  of  the  wing  and 
fuselage  must  be  included.  This  is  done  by  including  the  diffraction 
from  the  interior  wedge  illustrated  in  Fig.  38.  Note  that  once  the 
diffraction  point  is  known  only  the  wedge  angle  (which  is  defined 
by  "n"  in  the  diffraction  function)  is  changed  from  the  previous 
solution.  In  this  case  the  effective  source  is  found  by  using  the 
helical  path  solution  for  the  surface  rays  on  the  fuselage  except  now 
the  path  intersects  the  edge  formed  by  the  fuselage-wing  junction. 

The  tangent  direction  of  the  path  at  that  point  gives  the  incident 
field  direction  with  the  effective  source  '.ocation  for  the  diffracted 
field  calculation  being  shown  in  Fig.  38.  Mote  again  that  the  in¬ 
cident  field  values  are  actually  given  by  the  modal  solutions  of 
Section  II-B.  However,  in  this  case  the  field  is  incident  along  the 
surface  of  the  wedge;  thus,  the  incident  field  for  the  diffracted  field 
solution  from  this  edge  must  be  one  half  the  value  given  by  the  modal 
solution  as  discussed  in  Ref,  [493.  In  Fig.  39  the  radiation  pat¬ 
terns  for  the  three  infinitesimal  sources  are  shown  with  the  wings 
moved  up  such  that  =  45°  as  shown  in  Fig.  38.  The  validity  of  this 
solution  is  verified  by  the  measured  result  taken  for  the  monopole 
case  of  Fig.  39a. 

Even  though  the  previous  results  look  good  in  comparison  with 
the  measured  results,  several  improvements  in  the  solution  could  be 
made.  For  example,  it  was  assumed  that  the  fuselage  is  an  infinitely 
long  cylindrical  structure.  This  can  be  improved  by  including  the 
three-dimensional  effects  of  the  fuselage  as  introduced  in  the 
following  chapter.  Nevertheless,  based  on  data  taken  at  the  Naval 
Air  Development  Center  the  circular  shaped  fuselage  appears  to  be 
adequate  for  present  needs.  Other  possible  changes  that  might  be 
considered  are  tilting  the  wings,  mounting  the  antenna  off  the 
fuselage,  introducing  other  possible  scattering  structures  in  the 
analysis,  etc. 

C.  Elevation  Plane  Analysis 

This  section  is  basically  a  synopsis  of  the  material  presented 
in  Ref.  [48].  It  is  presented  here  to  illustrate  the  validity  of  the 
GTD  approach  for  treating  antennas  mounted  on  a  convex  two-dimensional 
body.  Note  that  this  material  is  an  essential  step  in  determining 
the  volumetric  pattern  of  an  antenna  mounted  on  a  three-dimensional 
surface  as  considered  in  the  next  chapter. 
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Fig.  39a— Roll  plane  pattern  of  monopole  (E  ). 


Based  on  our  original  aircraft  model  of  Fig.  28,  the  fuselage  is 
composed  of  cylinders,  cones,  and  spheres.  It  was  shown  by  Ryan  in 
Ref.  C20]  that  a  two-dimensional  diffraction  analysis  could  be 
adapted  to  solve  for  the  elevation  plane  pattern  for  an  antenna 
mounted  on  a  rocket  model  composed  of  cylinders  and  cones.  With  the 
antenna,  in  our  problem,  being  mounted  on  the  fuselage  near  the  top 
or  bottom  of  the  aircraft  one  can  to  a  good  approximation  just  con¬ 
sider  the  fuselage  effect  in  the  elevation  plane  in  that  the  other 
terms  have  normally  a  secondary  effect.  Thus,  Ryan's  approach  is 
applied  here  with  only  the  spherical  nose  cone  being  included  in 
his  solution  in  order  to  complete  our  analysis  of  the  originally 
specified  model. 

Some  of  the  results  of  the  present  study  are  illustrated  in 
Fig.  40.  In  Fig.  40a  Ryan's  calculated  result  is  compared  with  the 
the  measured  result  for  a  circumferential  slot  mounted  on  the 
fuselage  with  good  agreement  obtained  between  the  two  results. 


Fig.  40a— Elevation  plane  patterns. 
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Fig.  40b— Elevation  plane  patterns. 


In  Fig.  40b  a  monopole  is  placed  on  the  fuselage  with  the  analytic 
results  compared  with  the  isolated  dipole  pattern  showing  that  in 
the  lit  region  the  fuselage-mounted  antenna  radiates  much  like  the 
isolated  antenna.  However,  the  on-aircraft  antenna  pattern  deviates 
greatly  outside  the  lit  region.  Thus,  one  must  be  careful  when 
approximating  the  on-aircraft  antenna  performance  using  simply  an 
isolated  antenna  pattern.  Finally,  in  Fig.  40c  the  radiation  patterns 
of  monopoles  mounted  on  different  fuselages  are  considered  to  il¬ 
lustrate  the  effect  of  the  fuselage  shape  on  the  resulting  pattern. 
Note  that  with  the  spherical  nose  cone  the  ripple  is  greatly  reduced 
with  much  more  energy  radiated  in  the  forward  direction. 


Ill 


Fig.  40c— Elevation  plane  patterns. 


In  the  previous  fuselage  models  as  well  as  the  ones  to  follow, 
the  nose  radome  has  not  been  treated  as  a  special  case.  In  other 
words,  it  is  assumed  that  the  nose  radome  is  a  perfect  conductor  as 
is  the  rest  of  the  fuselage.  These  solutions  can  be  in  error  due 
to  this  approximation  especially  for  dielectric  radomes.  However, 
it  appears  that  metallic  radomes  will  find  widespread  use  in  the 
future;  in  which  case,  our  model  could  be  a  good  approximation  in 
the  radome  region  especially  outside  its  pass-band. 

The  simple  models  described  above  are  not  general  enough  to 
include  the  wide  variety  of  aircraft  fuselages  encountered  in 
practice.  An  aircraft  fuselage  is  predominantly  a  convex  body  which 
can  not  be  completely  described  by  simple  analytic  equations.  In 
practice,  an  aircraft  fuselage  is  often  specified  by  a  set  of  points. 
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Consequently  a  new  approach,  called  "Section  Matching  GTD  Method," 
is  developed  in  which  a  set  of  discrete  points  is  used  to  outline  the 
profile  of  the  fuselage.  In  this  way  any  convex  fuselage  shape  can  be 
included  in  our  general  solution  in  a  straight-forward  manner.  Mote 
that  this  solution  is  still  two-dimensional  but  in  many  cases  this  is 
a  good  approximation  in  computing  the  three-dimensional  elevation  plane 
pattern.  Nevertheless,  the  complete  three-dimensional  shape  and 
volumetric  pattern  will  be  treated  in  the  next  chapter  which  will 
remove  this  present  restriction. 

It  was  shown  in  Section  II-C  that  one  can  write  high  frequency 
asymptotic  expressions  for  the  far  zone  radiated  fields  of  an 
antenna  mounted  on  a  general  two-dimensional  convex  surface.  This 
solution  is  broken  up  into  solutions  for  the  three  regions  which 
ware  illustrated  in  Fig.  7.  7h'  forms  of  these  expressions  are 
given  in  general  terms  by  Eqs.  (22),  (24)  and  (25).  Note  that  the 
solution  in  the  lit  region  is  just  the  geometrical  optics  field  for 
the  isolated  antenna  and  is  not  affected  by  the  surface  geometry.  On 
the  other  hand,  the  solution  in  the  transition  region  is  dependent  on 
terms  of  the  form 


-j  k  ds 

e  •’  which  are  characteristic  of  the  convex  surface.  Note  that  in 
the  above  equations  pg(s)  is  the  radios  of  curvature  and  s  is  the 
arclangth  along  the  appropriate  geodesic  path.  If  the  surface  is 
now  approximated  by  a  set  of  points,  ore  can  evaluate  the  above 
integrals  numerically  using 


-j  [k  ds  l  ’'Asi  -jkas- 
e  J  2l  e  '  =  n  e 

i 


These  expressions  require  that  the  incremental  arclength  between  points 
(As-j)  and  the  radius  of  curvature  at  each  point  (pg)  be  found  based 
on  the  set  of  defining  points.  If  the  actual  surface  is  given  by 
f(x)  which  is  known  only  at  a  finite  number  of  points,  one  can 
make  a  polynomial  approximation  to  a  segment  of  the  surface  which 
is  given  by 
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The  geometry  for  this  problem  is  illustrated  in  Fig.  41.  Equating 
the  above  polynomial  solution  with  the  surface  equation  at  "n" 
consecutive  surface  defining  points  one  can  find  a  matrix  of  the 
form 


This  matrix  equation  can  then  be  inverted  to  solve  for  the  values 
of  the  c's.  Substituting  these  values  into  Eq.  (78)  one  obtains 
an  equation  which  locally  describes  the  surface  and  from  which  one 
can  determine  the  radius  of  curvature  and  incremental  arclength. 


Fig.  41 —Numerical  description  of  convex  surface. 


The  radius  of  curvature  can  be  defined  in  many  ways.  In  this 
case  it  is  given  by 


(79)  Pg  . 


which  is  applied  unless  dy/dx  approaches  infinity  in  which  case  let 
x  be  defined  in  terms  of  y  than 


The  values  of  the  various  derivatives  are  then  found  from  the 
polynomial  expression  of  Eq.  (78)  which  is  evaluated  at  the  mid  point 
(xn+1 )  of  the  approximating  section  by 
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giving 


surface  defining  point.  Note,  however,  that  this  requires  a  new 
polynomial  solution  (Eq.  (78) )  for  each  pg  calculation.  One  should 
realize  that  the  arclength  as  well  as  the  radius  of  curvature  must 
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be  determined  to  a  good  approximation  in  that  they  will  be  used  in 
the  deep  shadow  region  to  find  dpg/ds  and  d2pg/ds2.  Consequently, 
using  the  straight  line  distance  Between  adjacent  points  is  not  the 
best  method  to  find  the  incremental  arclength  based  on  our  limited 
comparisons.  These  comparisons  were  made  between  various  numerical 
solutions  and  the  known  solutions  for  s,  pg>  pg,  and  p'g  along  an 
elliptical  surface.  A  better  approximation  appears  tosbe  one  which 
takes  into  account  the  curved  nature  of  the  surface  between  points. 
One  such  approach  is  illustrated  in  Fig.  42  where 


(81) 


Asi 


=  Pgfoj)  +  Pg(Xj+l) 


A0i  = 


with 


W  +  °g(xi^):ls1n'1(p9(xi)  ♦  pg(«w)j 


‘  yi)<: 


^xi+l 


This  information  can  then  be  used  in  Eq.  (24)  to  determine 
the  values  of  the  radiated  field  in  the  transition  region.  Recall 
that  g*U)  and  g*U)  are  Fock  functions  which  were  defined  in 
Section  II-D  and  are  tabulated  in  the  literature[50].  These 
tabulated  values  can  be  applied  to  a  simple  interpolation  procedure 
that  can  be  employed  to  compute  the  complete  set  of  values,  which 
completes  the  approach  used  in  the  transition  region. 


In  the  deep  shadow  region  additional  information  is  needed  about 
the  surface.  For  the  diffraction  and  attenuation  constants  one 
needs  information  about  pg,  dpg/ds,  and  d2pg/ds2  at  each  point 
Cf(xi)3  along  the  surface.  The  values  of  tne  radius  of  curvature 
pq(x-j)  and  incremental  arclength  s-j  have  been  defined  by  Eqs.  (80)  and 
(81 ) .  These  values  can  now  be  used  to  compute  the  other  terms  using 
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Fig.  42~ Geometry  describing  the  arc  length  computation. 


and 
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which  are  numerical  averaging  solutions  specifically  designed  to 
compute  derivatives  based  on  discrete  data  as  presented  in  Ref.  [51]. 

The  only  undefined  terms  remaining  are  the  values  of  Aj(<L)  and 
^j(~9m)  which  are  tabulated  Airy  functions.  In  our  solutions  for 
the  GtD  fields  in  the  deep  shadow  region  only  two  modes  are  considered. 
The  values  of  these  terms  for  the  first  five  modes  are  given  in 
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Section  II-D.  Thus,  it  has  been  shown  that  the  GTD  solution  in  the 
deep  shadow  region  can  be  completely  defined  in  terms  of  our  general 
two-dimensional  convex  surface  which  is  only  defined  by  a  finite 
number  of  points. 

For  the  actual  computations,  the  tangent  values  at  each  of  the 
defining  points  were  determined  using  the  polynomial  approximation 
to  a  section  of  the  surface.  These  tangent  directions  in  turn 
specify  the  radiation  direction  of  the  fields  diffracted  from  that 
point.  Note  that  curved  surface  diffraction  is  a  local  surface 
phenomenon  as  shown  in  Ref.  [523.  Thus,  the  points  which  define 
the  surface  must  be  spaced  close  enough  together  such  that  the 
tangent  directions  of  adjacent  points  do  not  vary  too  rapidly. 
Otherwise,  the  pattern  can  only  be  determined  for  large  angular 
changes,  which  means  some  of  the  pattern  structure  can  be  lost  due 
to  the  large  pattern  change  between  data  points. 

It  was  found  for  the  cases  considered  in  this  study  that  a 
4th  order  polynomial  was  sufficient  to  locally  describe  the  surface. 
This  was  based  on  a  comparison  of  our  numerical  data  with  known 
results  for  various  convex  surfaces.  Note  that  the  values  of  pg, 
pg,  and  pg  could  then  be  computed  at  the  center  point  of  each  approxi¬ 
mated  section.  It  is  this  section  approximation  of  the  surface  that 
leads  to  the  term  Section  Matching  GTD  Solution.  Actual  numerical 
calculations  are  presented  in  Ref.  [483  to  show  the  validity  of  these 
various  approximations. 

In  order  to  verify  our  solution,  it  was  first  compared  with  the 
modal  solution  (Eqs.  (3)  and  (4))  for  an  antenna  mounted  on  a 
circular  cylinder.  These  solutions  are  compared  for  the  infinitesimal 
antennas  as  shown  in  Figs.  43  for  a  one  wavelength  radius  cylinder. 

In  each  case  there  is  very  good  agreement  between  the  two  results. 
These  results  do  tend  to  verify  this  approach.  However,  this 
fuselage  profile  is  circular  (pg  =  pg  =  0)  and  the  circular  case 
is  not  a  true  test  for  our  more  general  solution. 

This  section  matching  GTD  solution  now  is  extended  to  obtain 
radiation  patterns  in  the  elevation  plane  for  a  simulated  fuselage 
model  with  an  elliptical  profile.  The  calculated  radiation  patterns 
for  infinitesimal  antennas  such  as  circumferential  slot,  axial  slot 
and  infinitesimal  monopole  mounted  on  an  infinitely  long  elliptical 
cylinder  with  a  semi-major  axis  a  =  0.637a  and  semi-minor  axis 
b  =  0.390a  are  shown  in  Figs.  44  to  46.  The  comparison  between 
the  continuous  GTD  solution[373  and  the  section  matching  GTD 
solution,  again,  is  very  satisfactory.  In  addition  these  methods 
have  been  applied  to  more  general  fuselage  models  such  as  composite 
elliptical  cylinders.  The  radiation  patterns  in  the  elevation  plane 
for  a  circumferential  slot  mounted  at  various  locations  on  a  composite 
elliptical  cylinder  are  shown  in  Fig,  47.  The  results,  again,  compare 
very  favorably  with  the  continuous  GTD  solution[373.  The  elevation 
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plane  patterns  .'or  an  axial  slot  and  infinitesimal  monopole  are 
also  shown  in  Figs.  48  and  49,  respectively.  Based  on  these  results, 
the  pattern  in  the  lit  region  is  simply  the  direct  radiation  by  the 
antenna  itself.  On  the  other  hand,  the  pattern  in  the  shadow 
region  is  completely  dominated  by  the  fuselage  structure  through  the 
diffraction  mechanisms-  Thus,  the  complete  radiation  pattern  for  an 
antenna  mounted  on  a  fuselage  is  greatly  affected  by  the  fuselage  as 
one  should  expect.  The  location  at  which  the  antenna  is  mounted  is 
also  an  important  factor  in  determining  the  radiation  pattern.  The 
radiation  patterns  in  the  elevation  plane  for  a  smaller  composite 
elliptical  fuselage,  with  ai  =  1.5X  and  =  2. Ox  and  b  =  0.5X, 
are  presented  in  Fig.  50.  From  the  comparison  between  these  figures 
one  notices  that  there  are  fewer  backlobes  but  more  energy  radiated 
in  the  shadow  region  for  the  smaller  structure.  The  larger  structure 
has  more  backlobes  at  a  lower  average  energy  level.  This  is  simply 
due  to  the  size  of  the  fuselage  and  the  interaction  between  the  two 
surface  waves  propagating  in  opposite  directions  from  the  source;  in 
which  case,  the  two  ray  paths  play  a  most  important  role.  That  is, 
the  longer  the  ray  path,  the  more  the  surface  ray  attenuates. 

Again,  this  shows  the  significant  effect  which  the  structure  of  the 
aircraft  fuselage  plays  in  terms  of  the  antenna's  performance. 


D.  Azimuth  Plane  Analysis 

With  an  antenna  mounted  near  the  top  or  bottom  of  the  fuselage, 
it  is  quite  apparent  that  the  dominant  aircraft  structure  effect  in 
the  azimuth  plane  is  the  fuselage.  However,  in  this  case  one  is 
not  dealing  simply  with  the  cross-sectional  shape  of  the  fuselage 
as  in  the  roll -plane  or  the  profile  of  the  fuselage  as  in  the 
elevation  pKne.  In  fact,  in  order  to  deal  with  this  problem  to  a 
good  approximation,  a  more  complete  three-dimensional  study  must  be 
made  of  the  fuselage  as  is  considered  in  the  next  chapter.  On  the 
other  hand,  this  section  ir,  concerned  with  characteristic  features  of 
the  azimuth  pattern  ir.  terms  of  the  important  secondary  components. 

Two  such  contributors  considered  here  are  the  engine  and  wing 
scattering  effects. 

In  order  to  apply  our  previous  solutions  with  only  minor  modi¬ 
fications  it  is  again  assumed  that  the  antenna  is  mounted  on  an 
infinitely  long  circular  cylinder  (fuselage).  A*  in  the  roll  plane, 
the  modal  solutions  are  applied  to  determine  the  radiation  pattern  of 
the  fuselage  mounted  antenna.  Using  this  approach  one  can  again 
consider  an  arbitrary  antenna  by  integrating  the  equivalent  aperture 
currents  as  discussed  in  Section  II-B.  These  patterns  will  be  com¬ 
pared  with  patterns  including  additional  scattering  structures.  This 
should  allow  one  to  determine  if  these  secondary  scattering  structures 
are  important  in  the  azimuth  plane. 
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Fig.  44b--Elevat;on  plane  patterns  with  a  circumferential  slot 
mounted  on  an  elliptical  cylinder. 
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Fig.  47c— Elevation  plane  pattern  with  a  circumferential  slot 
mounted  on  a  composite  ellipse  model. 


Fig.  48a— Elevation  plane  pattern  with  an  axial  slot  mounted 
on  a  composite  ellipse  model. 
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Let  us  first  consider  the  engine  effect,  which  is  modelled  as 
a  finite  circular  cylinder.  Recall  that  in  the  roll  plane  the  engines 
had  little  effect  and  ware  neglected  in  the  more  general  three- 
dimensionfl  roll  plane  problem,  which  was  a  great  simplification  in 
the  analysis.  The  solution  of  this  problem  is  simply  an  extension  of 
the  finite  cylinder  scattering  problem  treated  in  Section  III-C.  The 
geometry  used  in  this  problem  is  illustrated  in  Fig.  51.  Note  that 
only  the  source  fields  incident  on  the  engine  have  been  modified  in 
this  situation  which  is  accomplished  using  the  techniques  developed 
in  the  roll  plane  analysis  of  Section  IV-B,  in  order  to  find  the 
effective  source  locations.  For  example,  the  source  field  incident 
on  the  rear  end  cap  of  the  engine  is  illustrated  in  Fig.  52.  The 
location  of  the  center  of  the  cap  is  given  by 

(xd  *  0,  yd  =  d,  zd  -z'+  a/2). 

With  this  location  known  one  can  apply  Eq.  (76)  to  find  the  ef¬ 
fective  source  location  on  the  fuselage  which  is  given  by 


ye  -  4/d 

X  (-ZM/2H* -*J 


where  <j>e 


With  the  effective  source  location 


known  one  can  apply  the  end  cap  scattering  effects  as  given  by 
Eqs.  (51),  (64).  The  source  field  value  incident  on  the  rear  end 
cap  is  given  by  the  modal  solutions  with  the  radiation  direction 
defined  by 


e  =  tan 


-1 


J(-*e)2  ♦  <i-y/  \ 

.  -  z‘-  ze  j 


4>  3  tan 


-1 
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TOP  VIEW 


Fig.  52— Geometry  of  source  field  incident  on  rear  end 


Similar  expressions  can  be  developed  for  the  effective  source 
location  and  field  value  for  the  front  end  cap. 

In  order  to  include  the  effect  of  the  cylindrical  section  one 
can  apply  the  geometry  illustrated  in  Fig.  53.  The  field  incident 
on  the  engine  is  dependent  on  the  desired  radiation  direction  (es, 
$s  =  90°)  and  determined  using  the  modal  solutions  with  the 
direction  given  by 


Again  the  results  of  Section  IV-B  can  be  applied  to  determine  the 
effective  source  position  which  in  turn  can  be  used  in  the  near  field 
cylinder  problem.  Using  Eq.  (77)  ore  finds  that 

xe  =  af  cos (90°-$) 

ye  =  af  sin(90%) 

z  =  af  cos (90°-$) tan  (9O°-0)(n/2  -$-$). 

with  the  actual  source  location  on  the  fuselage 
Zso  -  0).  Note  that  this  term  is  only  included 
of  the  cylindrical  section  as  shown  in  Fig.  53. 
in  the  cylinder  scattered  field  should  be  compensated  for  by  the 
edge  diffracted  field  contributions  from  the  ends  of  the  finite  engine. 
The  engine  scattered  field  is  then  added  to  the  directly  radiated  source 
field  using  the  superposition  principle  to  obtain  the  total  radiation 
pattern. 

Some  of  the  results  of  this  study  are  presented  in  Figs.  54.  Note 

that  in  each  case  the  engine  has  very  little  effect  on  the  total 

solution.  The  maximum  deviation,  for  the  cases  considered,  is  less 
than  2  dB  between  the  above  solution  and  the  solution  that  considers 
only  the  directly  radiated  term.  This  implies  that  one  can  to  a  good 
approximation  leave  the  engine  effect  out  of  the  solution  in  the 
azimuth  plane  provided  that  the  source  is  limited  to  the  region 
originally  assumed.  Finally,  measured  patterns  were  taken  for  the 
monopole  case  in  Fig.  54a  in  which  it  was  found  that  the  engine  played 
such  a  little  effect  that  it  was  difficult  to  observe  the  measured 
pattern  changes  with  and  without  the  engine  (finite  cylinder). 


defined  by  (af,^, 
over  the  finite  limits 
This  discontinuity 
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WITH  ENGINE 


The  wing  and  possibly  horizontal  stabilizer  effect  in  the  azimuth 
plane  must  also  be  considered.  To  analyze  these  structures  one  can 
simply  apply  the  solutions  presented  for  the  three-dimensional  roll 
plane  problem  of  Section  IV-B  with  only  the  desired  radiation  direction 
(9sj$s)  being  changed  to  those  angles  dictated  by  the  desired  azimuth 
pattern.  Again  these  results  are  compared  with  and  without  the  finite 
wing  effect  in  order  to  ascertain  their  significance  in  the  overall 
azimuth  plane  pattern. 

Some  of  the  results  of  this  study  are  presented  in  Figs.  55. 

Again  it  is  observed  in  each  case  that  the  wing  has  little  effect  on 
the  azimuth  plane  patterns.  In  fact,  what  little  effect  it  does 
have  occurs  in  the  small  sector  aft  of  the  aircraft  which  mi^ht  only 
be  significant  for  our  assumed  infinite  cylinder  model.  Thus,  the 
wing  can  be  neglected  to  a  good  approximation  in  the  azimuth  plane 
provided  that  the  assumed  geometry  and  source  location  are  adequate 
to  represent  the  actual  on-aircraft  antenna  problem. 

As  a  result  of  this  study,  it  is  apparent  that  the  dominant 
scatterer  in  the  azimuth  plane  is  the  fuselage.  This  effect  must  be 
taken  into  account  in  terms  of  a  three-dimensional  model  of  the 
fuselage  if  the  resulting  pattern  is  to  be  representative  of.  actual 
on-aircraft  antenna  patterns.  This  problem  is  considered  in  the 
next  chapter.  Nevertheless,  based  on  the  results  of  this  section 
one  can  to  a  good  approximation  compute  the  azimuth  plane  pattern 
simply  by  treating  the  antenna  mounted  on  a  reasonable  model  of 
the  fuselage  without  analyzing  the  other  secondary  scattering 
effects. 
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CHAPTER  V 


VOLUMETRIC  PATTERNS  OF  ANTENNAS  MOUNTED 
ON  CONVEX  SURFACES  OF  REVOLUTION 


As  stated  earlier,  for  a  source  mounted  on  the  fuselage  of  an 
arbitrary  aircraft  near  the  extreme  top  or  bottom  the  fuselage 
shape  has  the  dominant  effect  on  the  resulting  antenna  pattern.  The 
wings  and  other  flat  plate  structures  can  have  strong  effects  in 
certain  sectors  of  the  pattern  but  they  are  not  as  dominant  as 
the  fuselage  especially  when  the  complete  volumetric  pattern  is  con¬ 
sidered.  For  these  reasons  the  volumetric  pattern  of  an  antenna 
mounted  on  a  three-dimensional  isolated  fuselage  is  analyzed  by  the 
method  presented  in  this  chapter.  Nevertheless,  the  wings  and 
various  other  structures  could  be  added  to  this  solution  in  the 
future  as  was  done  in  tha  three-dimensional  roll  plane  problem. 

In  the  previous  chapter,  aircraft  models  were  considered  that 
would  resemble  a  wide  variety  of  aircraft  shapes  and  yet  could  also 
be  analyzed  with  reasonable  accuracy.  In  this  case,  it  is  quite 
obvious  that  the  three-dimensional  nature  of  the  fuselage  must  be 
modelled  if  one  is  to  adequately  determine  volumetric  patterns. 

In  the  elevation  plane  it  was  found  that  the  profile  of  the  fuselage 
had  to  be  accurately  represented.  This  resulted  in  the  development 
of  the  Section  Matching  GTD  solution  in  which  the  profile  was 
described  by  a  finite  number  of  points.  Whereas  in  the  roll  plane, 
the  circular  cross-section  was  found  to  be  adequate.  One  obvious 
extension  of  these  observations  is  to  consider  a  fuselage  which  is 
analytically  described  by  a  surface  of  revolution.  In  this  case 
the  profile  is  again  described  by  a  set  of  points  which  in  turn  are 
revolved  about  the  axis  of  the  fuselage.  In  this  way  the  important 
features  of  our  past  studies  in  the  three  principal  planes  are  in¬ 
corporated  in  this  new  three-dimensional  study.  Not  only  is  this 
shape  quite  versatile  but  it  can  also  be  analyzed  by  extending  the 
techniques  that  were  developed  previously.  This  will  become 
apparent  in  the  following  discussion. 

As  presented  in  Section  II-D  the  rays  which  propagate  outward 
from  a  source  travel  around  the  surface  along  geodesic  paths  while 
energy  is  continually  being  diffracted  in  the  tangent  direction  at 
each  point  along  the  path.  Thus,  the  first  obstacle  in  computing 
the  volumetric  pattern  of  an  antenna  mounted  on  a  three-dimensional 
convex  body  is  finding  a  numerical  technique  to  specify  the  geodesic 
paths.  One  such  solution  is  based  on  tensor  analysis  from  which  two 
differential  equations  are  given  in  Appendix  II  as 
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(84) 


R  sin  e(R  cos  e+  ^ 


sin  e) 


/dtf . 

(3lj  ' 


d2*  + 

37 


HP 

2(R  cos  e  +  ^  sin  e) 
R  sin  e 


de  d*  _  n 
HI  31"  0 


where  £  is  the  arclength  of  the  geodesic  path,  "i he  geometry  of  this 
problem  is  illustrated  in  Fig.  56.  Note  that  the  surface  is  defined 
by  a  set  of  points,  which  are  used  to  specify  R(e)  for  0  <_  9  <_  180°. 
For  the  cases  considered  here  R(e)  is  defined  every  %°  which  requires 
a  total  of  361  points  to  define  the  surface. 

Consider  a  point  on  the  surface  defined  by  e  and  $  which  is  a 
function  of  the  arclength  U).  Using  a  power  series  expansion  one 
finds  that 


o(0 


j  0(n)(0) 

n=0  n! 


t:11  ,  and 


*(i) 


£  *(m)(0) 

m=0  ml 


If  it  is  assumed  that  one  moves  a  very  short  distance  (a i)  along  a 
given  geodesic  path  then 


de 


e(i+Ai)  2:  eU)  + 


(85) 


A0  *  1  d6e 

*37 


A£  ,  and 


<f>  ( S.+A5.)  ^  (f>  ( £ )  +  — ^  A£  +  i  — —y 

di  i  c  dr 


A£.2  . 


Now  referring  back  to  thf  differential  equations  of  Eq.  (84)  one 
must  define  the  source  location  and  starting  direction  in  order  to 
completely  specify  a  geodesic  path.  However,  with  this  knowledge  one 
can  use  Eq.  (84)  to  find  the  second  derivatives  of  e  and  <j>  with  respect 
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to  £  at  the  point  (e,*)  as  given  by 


expressions.  Substituting  these  results  into  Eq.  (85),  the  values  of 
e(£+A£)  and  $(£+m)  are  located  approximately  which  in  turn  defines  a 
new  point  along  the  geodesic  path.  The  derivatives  of  e  and  $  at  the 
new  point  are  given  by 


Knowing  these  terms  one  can  again  use  Eq.  (86)  to  find  ar 

d2 . |  d£^  £+A£ 

— |  which  can  be  used  to  locate  a  third  point  along  the  given 

d£  *£+A£ 


geodesic  path.  By  continuing  this  process  one  can  completely  trace 
out  a  geodesic  ray  path  on  an  arbitrary  surface  of  revolution.  Some 
examples  of  geodesic  paths  on  surfaces  of  revolution  are  illustrated 
in  Figs.  57  and  58.  Note  that  the  geodesics  on  a  sphere  are  great 
circles  as  is  found  to  be  the  case  in  checking  our  formulation  and 
is  shown  in  Fig.  57. 


An  important  question  yet  to  answer  is  just  how  large  can  A£  be 
without  causing  our  solution  to  be  inaccurate.  This  is  a  difficult 
question  to  answer;  however,  some  knowledge  of  the  error  can  be 
found  using  Clairaut's  theoremC53]  which  states  that 

r  sin  a  =  constant 


for  a  given  geodesic  curve  on  a  surface  of  revolution.  Note  that  r 
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is  the  radius  of  the  parallel  passing  through  the  given  point  on 
„he  geodesic  path  and  a  is  the  solid  angle  between  the  geodesic 
direction  and  the  meridian  passing  through  that  point,  the  meridian 
being  the  profile  shaped  curve  that  is  revolved  about  the  fuselage 
axis.  Thus,  one  can  examine  the  error  by  the  variance  of  the  constant 
defined  above  as  he  progresses  along  a  given  geodesic  path.  For 
example,  on  a  prolate  spheroid  with  semi -major  and  semi -mi nor 
dimensions  given,  respectively,  by  4x  and  2x  and  with  Lt  =  0.Q05X 
it  was  found  that  Clairaut's  theorem  was  satisfied  to  within  3 
significant  figures  for  several  geodesic  paths  as  illustrated  in 
Table  IV. 


As  stated  earlier  R(e)  is  defined  every  for  0  <  e  <  180°; 

however,  as  observed  in  the  previous  equation,  not  onTy  Rle)  but 

dR  d^R  „ 

and  — p  must  be  known  over  the  complete  surface.  To  determine 

this  information  a  simple  4th  order  polynomial  was  found  that  passed 
through  5  consecutive  points  defining  the  surface.  The  derivatives 
of  R( e)  are  then  determined  at  the  center  point  of  each  section  using 
the  same  apprcach  presented  in  Section  IV-C.  The  values  of  these 
terms  between  the  surface  defining  points  are  found  using 
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where  e  is  an  arbitrary  position  parameter,  6]  is  the  closest  %°  value 
to  e  which  is  also  less  than  e  and  02  =  0-j  +  h°. 

A  second  solution  for  the  geodesic  paths  is  found  using  the 
calculus  of  variations.  Using  this  approach  one  forms  an  integral  of 
the  arclength  along  a  path  on  the  surface  which  can  be  solved  for  an 
extremum.  This  extremum  curve  is  simply  the  geodesic  path.  In  this 
case,  the  integral  of  the  arclength  is  given  by 
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which  has  an  extremum  given  53U 
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with  C  being  a  constant.  One  can  then  relate  C  to  the  initial  direction 
of  a  given  geodesic  curve  and  compare  this  result  with  the  previous 
solution. 


TABLE  IV 


CLAIRAUT1 S  CONSTANT  ON  A  SPHEROID 
USING  TENSOR  ANALYSIS  SOLUTION 


ii 

_ i 

o 

o 

o 

Y 

=  120° 

Clairaut 

Clairaut 

e 

4> 

Constant 

e 

$ 

Constant 

89.97513 

0.14106 

1.96962 

89.92838 

0.12405 

1.73205 

87.00310 

17.07253 

1.96961 

85.07056 

8.56457 

1.73199 

85.01104 

28.65405 

1.96961 

80.00671 

17.53570 

1.73193 

84.01189 

34.59238 

1.96960 

75.03504 

26.70357 

1.73188 

82.00377 

46.91133 

1.96960 

70.04465 

36.49352 

1.73184 

80.01470 

59.82879 

1.96960 

65.02831 

47.23272 

1.73181 

79.00000 

66.80015 

1.96960 

60.02371 

59.28953 

1.731/8 

78.01033 

73.92622 

1.96959 

55.03424 

73.37006 

1.73177 

76.01345 

89.65047 

1.96959 

51.51429 

85.24570 

1.73176 

y  =  140° 

Y 

=  160° 

Clairaut 

Clairaut 

0 

<*> 

Constant 

9 

4> 

Constant 

89.89027 

0.09207 

1.28557 

89,86540 

0.04899 

0.68404 

80.08774 

8.39953 

1.28533 

79.97284 

3.68367 

0.68382 

69.99167 

17.49768 

1.28513 

70.08763 

7.52630 

0.68362 

60.03110 

27.70494 

1.28496 

60.02801 

11.92162 

0.68247 

55.04909 

33.58741 

1.28491 

50.06772 

17.15761 

0.68336 

45.00158 

48.15512 

1.28483 

40.02862 

24.04913 

0.68325 

40.04995 

57.48482 

1.28479 

30.02665 

34.07472 

0.68300 

35.04354 

69.43960 

1.28473 

25.04869 

41.34835 

0.68273 

33.00055 

75.42781 

1.28470 

20.13786 

51.48320 

0.68226 

160 


As  illustrated  in  Fig.  59  the  values  of  the  position  parameters 
(e,4>)  for  the  given  curves  compare  extremely  well  for  the  two 
solutions.  Note  that  the  above  Integra’  solution  is  evaluated 
simply  using  Simpson's  rule.  In  addition,  it  is  shown  in  Fig.  60 
that  the  tangent  directions  (et»<H:)  of  the  various  geodesic  curves 
found  using  the  two  solutions  compare  very  favorably.  Recall  that 
the  radiation  direction  is  specified  by  the  tangent  to  the  geodesic 
path.  The  actual  solutions  used  to  compute  the  tangent  directions 
will  be  presented  later. 


Fig.  59— Comparison  of  geodesic  paths  using  tensor  analysis 
and  calculus  of  variations  solutions. 


Finally,  in  Table  V  the  values  of  Clairaut's  constants  are  shown 
for  various  values  of  C  using  the  above  integral  solution  for  the 
geodesic  path.  Again  this  constant  remains  unchanged  through  3 
significant  figures  for  the  cases  considered. 
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TABLE  V 


's 


C*_AIRAUT ' S  CONSTANT  ON  A  SPHEROID 
USING  CALCULUS  OF  VARIATION  SOLUTION 


II 

o 

°o 

Y  = 

120° 

Clairaut 

Clairaut 

0 

4 

Constant 

6 

A 

Constant 

89.97513 

0.00000 

1.96962 

89.92838 

0.00000 

1.73205 

87.00310 

16.98105 

1.96961 

85.07056 

8.67377 

1.73199 

85.01104 

28.38179 

1.96961 

80.00671 

17.48598 

1.73193 

84.01189 

34.1 5925 

1.96960 

75.03504 

26.62195 

1.731S8 

82.00377 

45.96819 

1.96960 

70.04465 

36.30131 

1.73184 

80.01470 

58.27634 

1 .86960 

65.02831 

46.80939 

1.73181 

79.00000 

64.69234 

1.96960 

60.02371 

58.55346 

1.73178 

78.01033 

71.33482 

1.96959 

55.03424 

72.18228 

1.73177 

76.01345 

85.51688 

1.96959 

51.51429 

83.44494 

1.73176 

y  =  140° 

_  _  _ 

Y  = 

160° 

Clairaut 

Clairaut 

0 

$ 

Constant 

e 

❖ 

Constant 

89.89027 

0.00000 

1.28557 

89.86540 

0.00000 

0.68404 

80.08774 

8.46014 

1.28533 

79.97284 

3.66851 

0.68382 

69.99167 

17.42221 

1.28513 

70.08763 

7.54001 

0.68362 

60.03110 

27.56524 

1.28496 

60.02801 

11.87875 

0.68347 

55.04909 

33.38832 

1.28491 

50.06772 

17.07832 

0.68336 

45.00158 

47.60200 

1.28483 

40.02862 

23.33207 

0.68325 

40.04995 

56.7692? 

1.28479 

30.02665 

33.60774 

0.68300 

35.04354 

68,30339 

1.28473 

25.04869 

40.67232 

0.68273 

33.00055 

73.91109 

1.28470 

20/3786 

50.59393 

0.68226 
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Using  the  geometry  illustrated  in  Fig.  61  a  set  of  test  curves 
is  defined  as  a  function  of  the  starting  direction  angle  (-,)  where 
0  _  ,  :  180.  Note  that  the  source  position  is  defined  by  [R(e0), 
o0,  io  '=  01.  The  vectors  e]  and  e 2  are  the  unit  tangent  vectors  to 
the  surface  which  are  defined  in  the  •:  =  0  plane  by 


The  starting  direction  of  any  aroitrary  ray  path  is  defined  by 
s  =  cos  y  e-j  +  sin  y  e ^  s 

where  these  terms  can  be  related  to  de/d;,  and  d*/d>.  by 


These  expressions  give  the  initial  conditions  used  to  solve  the 
differential  equations  of  Eq.  (84). 

Based  on  our  results,  solving  for  the  test  curves  with 
Ay  ~  Zk°  was  satisfactory  to  find  adequate  values  for  the  actual 
curves  used  in  the  pattern  calculation.  For  example,  the  test 
curves  are  shown  in  Fig.  62  for  a  source  mounted  at  the  top  of 
a  prolate  spheroid.  Then  any  desired  radiation  direction  can  be 
plotted  on  this  graph,  which  in  turn  can  be  used  through  an 
interpolation  procedure  to  determine  the  value  or  values  of  the 
starting  directions  (y)  from  the  source,  in  this  way  all  the 
dominant  ray  paths  can  be  traced  out  with  the  associated  field 
values  summed  in  the  far  field  to  give  the  complete  radiation 
pattern.  Note  that  the  test  curves  need  only  be  determined  for 
0  <_  y  <  180°,  since  for  a  surface  of  revolution  the  results  in 
the  otRer  half  space  are  simply  the  mirror  image. 
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Fig.  61— Starting  direction  of  geodesic  path  which  is 
defined  by  the  angle  (y). 


To  illustrate  the  validity  of  this  technique  the  resulting 
radiating  angles  are  presented  in  Table  VI  for  various  desired 
radiation  directions.  Note  that  as  the  errors  increase,  the 
resulting  amplitudes  of  the  field  values  are  decreasing  which  tends 
to  reduce  the  total  error  in  the  final  result.  This  results  from 
the  fact  that  the  error  increases  with  path  length  [i)  due  to  the 
numerical  computations  but  the  longer  the  path  the  more  the  field 
is  attenuated  due  to  e~al  in  the  GTD  field  expressions. 

Now  that  the  various  geodesic  paths  are  identified,  one  must 
determine  the  various  parameters  necessary  for  the  field  calculations 
using  the  GTD  approach  presented  in  Section  II-C.  First,  let  us 
consider  the  longitudinal  and  transverse  radius  of  curvature. 

In  order  £o  Refine  these  terms  consider  the  curvilinear  coordinates 
given  by  t,  n,  and  6  as  illustrated  in  Fig.  63.  Note  that  t  is 
the  unit  tangent  vector*  njs  £he  unit  normal  vector,  and  6  is  the 
unit  binormal  vector  (b  =  t  x  n).  The  radius  of  curvature  is  defined 
by 
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TABLE  VI 


RADIATION  DIRECTIONS  DESIRED  VERSUS  ACTUALLY  COMPUTED  TANGENT 
DIRECTION  OF  GIVEN  GEODESIC  CURVE  USING  INTERPOLATION  PROCEDURE 


6  .  .  .  =  90° 

desired 

edesired 

70° 

e 

comp. 

^desired 

6 

ycomp. 

e 

comp. 

^desired 

4comp. 

90.000° 

100° 

99.969° 

70.002° 

100° 

100.011° 

90.000 

no 

110. 053 

69  998 

no 

110.019 

90.000 

■ 

120.023 

70.002 

120 

120.016 

90.000 

129.992 

70.007 

130 

130.001 

90.000 

139.962 

69.990 

140 

139.960 

90.000 

150.046 

70.004 

150 

150.059 

90.000 

160.015 

69.997 

150 

159.945 

90.000 

169.985 

69.969 

— 

170 

169.925 

(88) 


1 


i 


where  t(z)  is  the  derivative  of  £  with  respect  to  the  arclength  (*.). 
The  unit  tangent  vector  to  the  geodesic  path  is  given  by 


(89)  t (A)  =  -d%^  =  77 


where  C(e,$)  is  the  position  vector  which  defines  the  geodesic  path. 
The  values  of  de/d?,  and  d$/d?.  are  defined  at  each  point  along  the 
path;  however,  dt/de  and  dS/d*.  must  be  found  using 

*f  A  \  A 

C(e,<+.)  =  R(e)  sin  e  cos  $  x  +  R(e)  sin  e  sin  ^  y  +  R(e)cos  e  z 


=  sin  e  cos  $  +  R  cos  9  cos  $ 


in  $jy  *  ^ 


R  cos  e  sin  <My  -5- |  cos  e  -  R  sin  e  Jz 


>)x  +  (h?  sin 

1  sin  ojz 


e  sin  6  + 


=  -R  sin  e  sin  <f>  x  +  R  sin  e  cos  <j>  y. 


Thus,  the  tangent  vector  is  defined  at  each  point  (e,<j>)  along  a  given 
geodesic  path  by  substituting  the  above  expressions  in  Eq.  (89). 

Since  the  diffracted  rays  radiate  in  the  tangent  direction,  this 
solution  defines  the  radiation  direction  at  every  point  along  a 
given  path. 

The  derivative  of  the  unit  tangent  vector  is  given  by 


d^C  ( deY  x  9  d2C  de  d*  ,  dC  d2e  , 
j  +  239TfHT  31  +  3e  ^  + 


d2C  de 


(TC  ( diY  *  df  d2<j> 

7W  • 


The  terms  found  in  the  above  equation  have  all  been  defined  earlier 
except  for 


9-*  o 

HR  \  A 

— 5-  =  (  sin  e  cos  $  +  2  cos  e  cos  c  -  R  sin  0  cos  $)x 

de^  Vde^  d6  J 


( d“*R  dR  \  ~ 

+  1 — ^  sin  e  sin  $  +  2  cos  a  sin  $  -  R  sin  e  sin  <jjy 

/  d^R  HR  \». 

+ 1  — £  cos  e  -  2  sin  e  -  R  cos  ejz  , 


— x-  =  -R  sin  e  cos  $  x  -  R  sin  9  sin  <>y  ,  and 

d«/ 


de  d<+> 


=  -sin 


^  sin  e  +  R  cos  ejx  +  cos  <^jj^  sin 


0  +  R  cos  8 jy 


All  of  these  terms  are  then  determined  once  the  geodesic  path  is 
traced  out  to  the  point  u(e,$)  on  the  surface.  The  longitudinal  radius 
of  curvature  is  then  simply  obtained  by  substituting  these  quantities 
into  Eq.  (88). 

In  order  to  solve  for  the  transverse  radius  of  curvature,  one 
must  find  similar  relations  for  the  curve  on  the  surface  of  revolution 
which  is  orthogonal  to  our  given  geodesic  path  at  the  point  t(e,$). 

In  so  doing  the  geodesic  path  orthogonal  to  our  original  path  must 
first  be  found.  This  will  allow  us  to  use  the  previously  derived 
equations  which  in  this  case  give  us  the  transverse  radius  of 
curvature.  The  unit  tangent  vectors  at  any  point  on  the  surface  are 
given  by 


W 


sin  0  +  R  cos  0  cos  <p  x  + 


^^sin  8  +  R  cos  ojsin  <py  +  cos  e  -  R  sin  ejz  ,  and 


=  -sin  $  x  +  cos  <py 
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The  unit  tangent  vector  to  the  original  geodesic  path  (t)  and  the 
transverse  path  (£)  can  be  written,  respectively,  as 


t  =  t-j  e-j  +  tg  ^2  ,  and 

A  A  \ 

b  =  b1  e-j  +  b2  e2 


Since  t  •  b  =  0,  one  finds  that 

*  „  „  dC. 

b  =  -t2  e!  +  t]  e2  =  dT~ 


where  £{.(e,$)  ^e  position  vector  of  the  transverse  geodesic  path. 
Substituting  Eq,  (91)  into  the  above  expression  it  is  seen  that 


However,  dCt/ds>  is  defined  by  Eq.  (89)  which  when  substituted  into 
the  above  equation  gives 


Note  that,  the  values  of  ti  and  t2  are  defined  at  each  point  along 
the  original  geodesic  path  by  Eq.  (89).  With  these  quantities  known 
one  can  f i nd  d2e/d£2jtrans  and  d24>/de2[trans  using  the  geodesic 
solutions  of  Eq.  (86).  Knowing  this  information  the  transverse 
radius  of  curvature  is  simply  found  as  before  using  Eq.  (88)  with 
the  transverse  path  terms  being  used  in  this  case. 
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The  accuracy  of  this  approach  is  illustrated  by  a  comparison  with 
various  analytic  solutions  for  the  longitudinal  and  transverse 
radii  of  curvature.  Our  approximate  values  were  compared  with  the 
known  values  for  a  sphere  and  found  to  be  in  agreement  through  5 
significant  figures.  Our  values  are  compared  with  the  analytic 
solutions  on  a  prolate  spheroid  along  a  meridian  geodesic  path  in 
Table  VII.  Good  agreement  is  obtained  between  the  two  solutions 
for  the  cases  considered.  Note  that  the  derivatives  of  the  rad’ us 
of  curvature  are  not  included  in  this  study  to  date  in  that  the 
numerical  complexity  would  be  too  great  to  warrant  its  insertion 
at  this  time.  However,  a  method  of  computing  sucn  quantities  has 
been  presented  in  Section  IV-C  so  that  it  could  be  introduced 
if  desired  at  a  later  date. 


TABLE  VII 

COMPARISON  OF  ACTUAL  AND  COMPUTED  VALUES  OF  LONGITUDINAL  AND 
TRANSVERSE  RADII  OF  CURVATURE 


Prolate  Spheroid  (4x  x  2x) 

0 

pt  actual 

pt  comp. 

pg  actual 

pg  comp. 

n.o 

1.00 

1.00 

1.00 

1.00 

10.0 

1.15 

1.15 

1.54 

1.54 

20.0 

1.43 

1.43 

2.91 

2.93 

30.0 

1.65 

1.66 

4.47 

4.58 

40.0 

1.79 

1.80 

5.76 

5.85 

50.0 

1.88 

1.89 

6.69 

6.78 

60.0 

1.94 

1.94 

7.32 

7.36 

70.0 

1.98 

1.98 

7.71 

7.75 

80.0 

1.99 

2.00 

7.93 

7.97 

90.0 

2.00 

2.00 

8.00 

8.02 

— 

Another  quantity  necessary  in  our  calculations  is  the  spread 
factor  (J  dipo/dip) ,  which  is  representative  of  the  amount  of  separation 
between  adjacent  geodesic  paths  as  they  propagate  around  the  surface. 
This  type  of  separation  is  illustrated  in  Fig.  64.  This  parameter 
is  obtained  numerically  by  considering  two  adjacent  rays  initially 
separated  by  d  0,  then  determine  d  which  is  the  included  angle  between 
the  adjacent  tangent  vectors  at  the  point  of  diffraction.  Thus,  if 
tangent  directions  at  the  point  of  diffraction  are  given  by  L  and 
tg,  than 


Fig.  64— Illustration  of  the  spread  factor  (N<  dipQ/ dip)  terms. 


A  study  was  made  to  determine  just  how  large  d^o  could  be  made 
without  greatly  distorting  the  value  of'/dijWdi  at  the  diffraction 
point.  It  was  found  that  dip0  could  be  as  large  as  2h°  which  meant 
that  the  values  of dip0/dip  could  be  computed  and  stored  along  with 
the  test  curves.  Consequently,  one  need  only  specify  the  radiation 
direction  which  in  turn  specifies  the  initial  directions  (y)  of  the 
geodesic  paths  and  the  associated  spread  factors  (v|  dip0/dg>) .  These 
data  are  stored  at  the  beginning  of  the  program  in  that  they  are 
characteristic  of  the  body  and  source  location  being  considered. 

The  far  field  patterns  can  now  be  computed  using  the  previously 
defined  terms  which  provide  information  about  the  convex  surface  of 
revolution  under  consideration.  Since  the  antenna  is  assumed  to 
be  mounted  on  the  fuselage  of  the  aircraft,  one  needs  only  consider 
two  general  types  of  antennas.  These  being  the  monopole  type  which 
has  a  normal  component  electric  current  with  respect  to  the  surface 
and  a  sloe  type  antenna  which  has  a  tangential  component  magnetic 
current.  It  i.;  assumed  that  the  fields  launched  by  infinitesimal 
antennas  follow  the  solutions  specified  for  the  two-dimensional 
problem  of  Section  II-D.  Note  that  in  this  case  the  three-dimensional 
geometry  is  introduced  in  terms  of  the  geodesic  paths  and  associated 


longitudinal  and  transverse  radii  of  curvature  which  appear  in  the 
diffraction  and  attenuation  coefficients.  It  is  further  assumed  that 
the  normal  and  tangential  component  boundary  layer  fields  propagate 
around  the  surface  independently.  This  approximation  might  be 
rather  poor  if  the  torsion  along  the  geodesic  curve  varies  too  greatly; 
however,  this  point  is  not  well  understood  at  this  time.  In  any  event, 
the  results  presented  here-  will  be  compared  with  actual  measurements 
to  illustrate  the  validity  of  these  assumptions.  On  the  other  hand, 
torsion  is  a  surface-curve  relation  which  can  be  computed  using 
techniques  similar  to  the  previously  derived  relations. 

At  this  time  only  a  convex  body  is  considered  in  order  that  the 
radiation  direction  can  be  simply  defined  by  the  geodesic  tangent 
direction  as  was  indicated  earlier.  However,  the  study  of  a  concave 
body  is  an  important  topic  worthy  of  future  consideration. 

Let  us  first  investigate  the  monopole  antenna  type  whose  surface 
rays  propagate  outward  in  all  directions  from  the  source  with  a 
normal  component  E-field  (or  follow  the  hard  boundary  condition). 

The  following  is  a  summary  of  the  solutions  in  the  three  different 
regions  for  the  monopole  case  using  the  geometry  illustrated  in  Fig. 

65: 


Fig.  65— Geometry  of  monopole  problem. 
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Lit  Region 


(92)  E(o,i)  =  -  sin  om  om  F(source) 


Transition  Region 
a)  Lit  side 


(93)  E  =  h^in  em)%  ejk*g 


b)  Shadow  side 


{'  ?ferqr)1/3MJ  "  •F(ia"!|ent) 


(94)  E  =  ni 


di> 


o 


e  g’ 


?(w7r)1/3^i]j‘  F(tanger,t) 


Deep  Shadow 


(95) 


where 


E  =  I  n.  eJ  F,  (tangent) 
j  J  J  J 


E h  -t—  D  ne  "  1  1 

J3T  m=o  m  m  i 


Note  that  the  summation  over  "j'1  in  the  above  expression  indicates 
that  several  terms  can  contribute  in  the  deep  shadow.  An  example 
of  this  situation  is  illustrated  in  Fig,  66  where  four  rays  con¬ 
tribute  to  the  far  field  pattern.  The  only  term  in  the  above 
expressions  that  is  yet  to  be  defined  is  the  launch  coefficient 
which  is  given  by[353 


(96) 


in 


at  source  location. 


174 


V 


Fig.  66—The  four  dominant  GTD  terms  that  radiate  at 
(e=90°,  4>=1 45° ) . 


A 

Note  that  n  is  the  unit  normal  at  the  tangent  (or  diffraction)  point 
of  a  given  component  of  the  far  field  pattern,  and  F  (  )  is  a 
phase  factor  used  to  reference  the  phase  to  the  center  of  the  surface 
of  revolution. 

In  order  to  obtain  measured  patterns  off  the  principal  planes 
using  a  conventional  pattern  range,  which  only  revolves  about  a 
vertical  axis,  the  body  is  first  rotated  by  an  angle  (c)  as  shown  in 
fig.  67*  Thus,  as  the  body  turns  about  the  z'-axis,  one  obtains  the 
e'  and  4. *  components  of  the  radiated  field.  In  Fig.  68  the 
component  of  the  field  is  illustrated  with  £  =  0°  for  a  short 
monopole  mounted  on  a  4a  x  2x  prolate  spheroid.  For  this  principal 
plane  pattern  the  compari son A between  the  measured  and  calculated 
results  is  quite  good.  The  a1  and  4'  components  of  the  radiated 
field  for  the  same  body  are  illustrated  in  Figs.  69  and  70  for 
K  =  200  and  40°,  respectively.  In  these  off-principal  planes  casas 
the  agreement  between  the  results  is  again  quite  good. 
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Fig.  67— Rotation  of  convex  surface  in  order  to  obtain 
off-principal  plane  patterns. 


Even  though  the  above  problems  considered  seem  to  be  rather 
special  cases  of  the  more  general  problem  originally  treated  at  the 
outset,  the  results  are  encouraging  especially  for  the  off-principal 
plane  cases.  They  do  tend  to  verify  this  approach  and  the  associated 
approximations  at  least  for  the  monopole  cases  considered. 

The  next  case  to  be  considered  is  that  of  an  arbitrary  oriented 
slot  mounted  on  a  convex  surface  of  revolution.  It  was  shown  in 
Section  II-D  that  for  the  two-dimensional  problem  with  the  slot  mounted 
parallel  (orthogonal)  to  the  cylinder  axis  that  the  slot  radiated 
according  to  the  hard  (soft)  boundary  condition  in  the  principal 
plane.  If  an  arbitrary  orientation  of  the  slot  were  considered 
then  one  could  assume  that  the  pattern  is  given  in  the  principal 
plane  by 


(97) 


E  =  Esoft  sin  6  +  Ehard  cos  6 
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MEASURED 


MEASURED 


MEASURED 


MEASURED 


using  the  geometry  illustrated  in  Fig.  71.  Thus,  for  a  slot  antenna 
mounted  arbitrarily  on  a  three-dimensional  surface  with  a  volumetric 
pattern  desired  one  can  extend  this  approach  by  considering  that  the 
boundary  layer  field  has  a  tangential  component  given  by 


Etang  =  Esoft  C0S^-S>  =  Esoft  cos(-,-e)b 


and  a  normal  component  given  by 


^normal  '  Ehard  s1n!v'e)  '  Ehard  =1ntv-e)n 


where  y  is  the  starting  direction  of  the  geodesic  path  defined 
earlier.  Note  that  the  above  solutions  agree  with  Eq.  (97)  in  the 
principal  plane  (y=  90°).  It  is  then  assumed  that  these  two 
components  propagate  around  the  surface  independent  of  one 
another  according  to  the  formulas  specified  in  Section  II-D. 


FINITE  SLOT 


Fig.  71— Finite  slot  geometry. 


A 
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Using  these  results  one  can  define  the  slot  far  field  pattern  in  teras 
of  the  three  regional  solutions  by: 


Lit  region 


(S8)  E  =[e(cos  $  sin  e)  +  $(sin  g  cos  e  sin  >,  -  sin  e  cos  g)] 


F(source) 


Transition  region 
a)  Lit  side 


(99)  i =jn [*  jkv(:  j/3 

v.  S  1 

b[j/2 

V.  \2p  (£-j)/  -'V  9  /at  source 


cos(Y+g) 


F(tangent) 


bj  Shadow  side 


*  "Is  [^J5e'dkv((^-  )1/3^)s1n(v-6)] 


2V*i> 


s  e'iU9‘ 


(&>  Wfe) 

y  i  ' 


1/3  "I 

cos(y-b)] 
at  source 


F(tangent) 


Deep  Shadow  Region 


(101)  E  *  l  [0^  eJ  sin(Y-g)  +  bj  Ej  cos(Y>=g)UF(tangent) 
J 


l.  I  do  1  .  . 

h  J~2-  l  lie  *  1  ' 

x|d<  0=0  0  m  i 


K>  1  s  ,  -vJ  '-•)«, 

"Jar  1  i*  "* 


m=0  m  m  i 


Again  the  summation  over  "j"  in  the  above  expression  indicates  that 
several  terms  can  contribute  in  the  deep  shadow  region  as  shown 
previously  in  Fig.  66.  The  launch  coefficient  for  the  hard  boundary 
condition  is  given  by  Eq.  (96);  whereas,  the  soft  boundary  condition 
launch  coefficient  is  given  by[35] 


(102) 


'm 


Ai(-V- 


at  source  location  . 


Note  that  n  and  b  are,  respectively,  the  normal  and  bi normal  to  the 
geodesic  curve  at  the  point  of  diffraction  and  F(  )  is  simply  a 
phase  factor  to  i  efer  the  phase  to  the  center  of  the  surface  of 
revolution. 


In  order  to  measure  the  pattern  of  a  slot  antenna,  a  finite 
length  slot  is  considered.  Actually  for  our  measured  patterns  an 
open-ended  waveguide  is  used  except  the  width  is  cut  down  to  ap¬ 
proximately  one-tenth  of  the  wavelength.  For  this  antenna  one 
can  compute  the  pattern  by  numerically  integrating  the  patterns  from 
three  infinitesimal  slot  antennas  which  approximate  the  aperture 
fields  as  shown  in  Fig.  71.  This  approach  is  described  in  Ref.  [22] 
and  successfully  applied  in  Ref.  [47].  Using  this  solution  the 
computed  principal  plane  pattern  [c=0°  in  Fig.  67]is  compared  with 
the  measured  result  as  shown  in  Fig.  72  for  a  circumferential  slot 
(e=90O).  Note  that  these  results  are  again  taken  on  a  prolate 
spheroid.  The  calculated  and  measured  off-principal  plane  patterns 
for  c=20°  and  40°  are  shown,  respectively,  in  Figs.  73  and  74.  In 
each  case,  good  agreement  is  obtained  between  the  measured  and  cal¬ 
culated  results. 


Again  the  slot  cases  treated  are  rather  idealized  as  compared 
to  the  more  general  problem  originally  outlined.  However,  off- 
principal  plane  radiation  patterns  from  antennas  mounted  on  three- 
dimensional  surface  which  are  large  in  terms  of  the  wavelength  have 
not  been  analytically  computed  with  much  success  to  date.  Thus, 
one  must  begin  by  treating  specific  cases  which  are  easily  modelled 
and  computed.  In  this  way  one  is  able  to  ascertain  the  validity  of 
his  solution.  The  above  results,  then,  indicate  that  the  high 
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frequency  solutions  of  Kpuyoumjian[35]  can  be  applied  at  least  to 
this  idealized  problem.  One  must,  now,  extend  these  solutions  to 
determine  for  what  type  fuselage  "mdels  they  tend  to  break-down  if 
in  fact  they  do. 


CHAPTER  VI 
CONCLUSIONS 


High  frequency  solutions  for  the  radiation  patten, «  of  general- 
type  on-aircraft  antennas  mounted  on  the  fuselage  is  the  object  of 
this  research.  It  has  been  shown  previously  that  the  modal  solutions 
are  quite  adequate  fur  computing  these  patterns  provided  the  structure 
is  not  illuminated  too  strongly.  However,  modern  applications 
require  antennas  which  can  strongly  interact  with  the  aircraft 
structure.  Thus,  improved  analytic  models  must  be  applied  in  order 
to  develop  useful  theoretical  solutions.  These  solutions  can  then  be 
used  to  compute  the  antenna  design  and  location  for  a  given  appli¬ 
cation. 

The  basic  approach  applied  here  is  to  add  to  the  modal  solutions 
by  including  additional  structural  components  in  the  analytic  model. 
This  is  accomplished  using  a  ray  optics  technique  in  which  the 
modal  solutions  are  written  in  a  ray  form.  Using  this  scheme,  one 
is  able  to  cn'-ider  isolated  structural  components,  which  have  been 
removed  from  the  aircraft,  and  then  add  it  to  the  model  simply  by 
adjusting  the  field  incident  on  the  structural  scatterer. 

Two  near  field  scattering  problems  are  initially  studied  in 
Chapter  III.  These  are  a  finite  flat  plate,  which  is  used  tu 
simulate  the  wing,  and  a  finite  cylinder,  which  is  used  to  ap¬ 
proximate  the  .jet  engine^.  The  flat  plate  scattering  is  solved 
using  three-dimensional  wedge  diffraction  techniques.  Various 
measured  results  are  presented  to  verify  our  theoretical  solutions. 

The  finite  cylinder  scattering  is  solved  using  physical  optics,  the 
Geometrical  Theory  of  Diffraction,  and  modal  solutions.  Again 
measured  results  are  presented  to  illustrate  the  validity  of  this 
solution. 

These  near  field  scattering  solutions  are  then  added  to  analytic 
aircraft  models  to  compute  the  principal  plane  patterns  in  Chapter  IV. 
In  the  roll  plane,  it  is  found  that  the  jet  engine  has  little  effect 
for  a  two-dimensional  model  and  is,  thus,  neglected  in  the  more 
general  three-dimensional  study.  The  resulting  roll  plane  model  is 
basically  an  infinitely  long  circular  cylinder  to  which  wings  and 
horizontal  stabilizers  can  be  attached.  Note  that  the  wings  and 
horizontal  stabilizers  are  approximated  by  arbitrarily  many  sided 
flat  plates.  The  elevation  plane  model  approximates  the  profile 
of  the  fuselage  by  a  set  of  discrete  points  with  the  resulting 
surface  being  convex.  The  fuselage  is  found  to  be  the  dominant 
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scatterer  in  the  elevation  plane.  Finally,  the  azimuth  plane  model 
consists  of  an  infinitely  long  circular  cylinder  to  which  wings 
and  engines  can  be  added.  However,  it  is  shown  the  wings  and 
engines  have  such  a  small  effect  that  they  can  be  neglected  to  a 
good  approximation  in  the  azimuth  plane.  In  each  case,  measured 
results  were  taken  in  order  to  verify  the  various  solutions. 

These  solutions  provide  a  high-speed  analytic  tool  for  determining 
the  design  and  location  of  antennas  based  on  their  principal  plane 
performance.  For  example,  the  programs  delivered  to  the  Nava-  Air 
Development  Center  typically  run  a  pattern  in  30  seconds  or  less 
on  their  CDC  6600  digital  computer. 

An  additional  feature  of  these  solutions,  which  can  be  very 
important  in  certain  critical  cases,  is  that  one  can  trace  out 
the  dominant  pattern  terms  as  they  are  scattered  by  the  various 
aircraft  structure*.  In  tnls  way  one  can  quickly  ascertain  which 
structural  scatterers  are  distorting  the  pattern  in  a  critical 
region.  This  g^ves  him  the  option  of  taking  corrective  action  by 
changing  the  structure  or  by  properly  placing  absorber.  So  these 
solutions  not  only  provide  fast  pattern  computations,  but  they, 
also,  provide  the  antenna  designer  with  a  means  of  analytically 
considering  several  alternatives  to  improve  the  antenna's  per¬ 
formance. 

The  final  problem  considered  in  Chapter  V  is  the  volumetric 
pattern  analysis  of  an  antenna  mounted  on  a  convex  surface  of 
revolution.  This  model  describes  the  fuselage  profile  by  a  set 
of  points  which  are  revolved  about  the  fuselage  axis  in  order  to 
analytically  approximate  a  complete  surface.  The  wings  and  jet 
engines  are  not,  presently,  considered  in  this  study  in  that  the 
fuselage  scattered  fields  tend  to  dominate  the  overall  pattern  as 
determined  by  the  principal  plane  studies.  Nevertneless,  these 
structural  components  can  be  considered  in  the  future  as  was  done 
in  the  principal  planes.  The  volumetric  patterns  of  monopoles  and 
slots  using  the  asymptotic  high  frequency  solutions  of  KouyoumjianC35I] 
are  numerically  derived.  The  geodesic  paths  are  found  and  verified 
as  well  as  the  radii  of  curvature,  etc.  The  principal  and  off- 
principal  plane  patterns  are  computed  and  compared  with  measurements 
made  on  a  prolate  spheroid.  In  each  case,  good  agreement  is  obtained 
which  tends  to  verify  the  numerical  solutions. 
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Now  using  the  cylindrical  separable  fieldsQ54]  from  Table  VIII 
one  finds  that 
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The  form  of  the  scattered  field  can  now  be  postulated  giving 


„  jkz  cos  0.  «  fo\, 

E*  -  S  e  1  I  a  H<2)  (kP  sin  8.)e‘w 

*  n=-» 


where  the  an's  are  unknowns.  Using  the  boundary  conditions  at  the 
^urfice  of  the  cylinder,  at  which  the  tangential  component  of  the 
E-field  must  vanish,  one  finds  that 

E*  (p=a)  =  j  eJkZ  C°S  0i  *  l  [-  jn  j;(ka  sin  9.)  +  a 
jl  *  n=_c  n  in 

* 


,(2)  * 
n 

which  implies 


H*;2Mka  sin  e.)]  e"jn*  0 


0‘Cka  sin  e,) 

’n"1 


'  (ka  sin  e.) 


The  total  field  is  now  specified  by  the  above  relations  in 
conjunction  with  Table  VIII  and  is  given  by 
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jkz  cos  Q. 


J'(ka  sin  0.) 


E*  s  r  cin  ft —  I  2njn+1  ■Uto  sin  0,-)  -nrrr'  ~v 

ip  kp  sin  e-  p=1  (_nv^  r  u(2)' 


1  (ka  sin  e.) 


H^(kp  sin  6.) | sin  n<f> 


-I 


r*  t  » 

Ex*  =  "J  6 


jkz  cos  e. 


.  r  r  J  ‘(ka 

1  ^  "n  j"  Jn(kp  sin  9i}  •  nW 

n~0  n  L n  1  Hu; 


J^(ka  sin  0.) 


n  (ka  sin  0.) 


(104) 


1 2) 1 

H^1 ’  (kp  sin  0n.)  |  cos  n<f> 


E*  =  0 

J.Z 


+  j  cos  0.  e 

Hz  = - 1— 

xp  ZT 


jkz  cos 


«  r  J' (ka  sine.) 

I  ejnb'(kp  sin  ei)  "  17P - 

n=0  n  Ln  1  h[Z}  (ka  sin©.. 


n  (ka  sine.) 


(kp  sin  0.)j  cos  n<{. 

jkz  cos  0. 

+  cos  0.  e  00  nxiT  J ’ (ka  sin  0.) 

— Ro~z — sin  e -  1 2n  j  rMkp  sln  8i^  — nrr - 

A*  kp  ZQ  sin  ei  n=1  L  n  1  nj2)(ka  sin  0.) 


H^(kp  sin  e.)|  sin  n<j> 


->] 


sin  0.  jkz  cos  0. 


J^ka  sin  0.) 


?  £  enJ  lVkpsin  ei)  "  ~72T 

n=0n  Ln  i  HiZ)  (ka  sin  0,) 


;in  e^jji 


H^(kp  sin  ejjcos  n<f>  , 


where 


-of 
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i 


1 


e 


n 


n=0 

( 

otherwise. 


Note  that  these  solutions  give  the  total  field  at  every  point  in 
space  for  a  plane  wave  with  a  perpendicular  component  E-field  in¬ 
cident  on  an  infinitely  long  cylinder  of  radius  (a). 


The  total  field  for  the  parallel  component  E-field  incident  is 
found  in  the  same  way.  In  this  case,  the  incident  field  is  given 
by 

+-•  .  *  jk(z  cos  e.  +  x  sin  e.) 

E(j  =  (x  cos  -  z  sin  e.)  e  1  1  . 


Using  Eq.  (103)  it  is  easily  shown  that 


E  =  -sin  e.  e 
ii  z  i 


jkz  cos  e.  » 


'  I  jn  Jn(kp  sin  e.)e 


*jn<f> 


n=-c° 


Again  the  form  of  the  scattered  field  can  be  postulated  as 


E  ,  =  -sin  0.  e 

IlZ  1 


jkz  cos  0.  ~ 


1  l  *„  j"  H<2>(kp  sin  epe-J"*. 


n=-« 


The  total  tangential  E-field  must  vanish  on  the  perfectly 
conducting  cylindrical  surface  in  which  case 


.  jkz  cos  e.  »  n 

E.f-»(p=a)  =  -sin  0,-  e  1  r  -n' 

n=oo 


'HZ' 


i 


I  jnDn(ka  sin  0^  + 


an  H^2)(ka  sin  0.)]e"Jn<f  =  0 


which  implies 


& 


Using  these  results  along  with  the  separable  fields  of  Table  VIII  < 
finds  that 


l  jkz  cos  e.  « 

E  =  -0  cos  o,  e  1  l  c  j"Cj;(kp  sin  e  )  - 

n=0  "  1 

Jn(ka  sin  e.)  (2), 

'  i^('ka  „„  e.)  "»  (kp  sin  cos  no 

jkz  cos  e. 
t  cos  ei  e  1 

Ene  kp  sin  e. -  I  2n  J  Wn(kPsin  e.)  - 

i  n=l  n  1 

J  (ka  sin  e  ) 

105)  i™(ka  sin  ,  )  ""  (k°  Sln  6i)]  Sln 


Vz  =  -sin  0i  e 


jkz  cos  e. 


I  enjn[Un(kpsin  6.)  - 

n=n  "  "  I 


Jn(ka  sin  0.)  .  . 

-T^r—  : - -Hj  '(kp  sin  e.)j  cos  n$ 

(ka  sin  e^)  "  1 


jkz  cos  0.. 

Ht  =  .  e _ y  9n  .n+lri  ,  Jn(ka  sin  e.) 

"p  kpZ  sin  0.  2n  J  W  (kp  sm  0.)  — pj-x - 

0  1  n=1  H^;(ka  sin  0.) 


H<2!kp  sin  9^)3  sin  n<f> 


4.  _•  jkz  cos  e.  * 

Hz=i~e  1  V  e  inrj» 

“*  ^  Jo  I3  L " 


V  <nrii/t  •  x  Jn(ka  sin  0i> 

1  EnJ  CJ'(kP  sin  0.)  -  JLr - 1 — 

=°  "  1  H^tka  sin  e.) 


(2) » 

(kp  sin  e^jcos  n* 


E 


Note  that  these  solutions  give  the  total  fields  at  any  point  in 
space  for  a  parallel  component  plane  wave  field  incident. 

The  solution  for  an  arbitrary  plane  wave  field  incident  can  be 
easily  found  by  decomposing  the  field  incident  into  its  perpendicular 
and  parallel  components  multiplying  these  values  by  the  above 
solutions  and  sunning  the  terms. 

It  was  assumed  for  convenience  in  the  above  modal  solutions  that 
the  plane  wave  was  incident  at  $i  =  0;  hov/ever,  this  restriction  is 
simply  removed  by  replacing  $  in  Eqs.  (104)  and  (105)  by  $-$.  where 
$i  is  the  general  incident  $  direction. 
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TABLE  VIII 

TE  AND  TM  FIELDS  SEPARABLE  IN  THE  CYLINDRICAL  COORDINATE  SYSTEM 

The  harmonic  electromagnetic  fields  listed  below  satisfy  Harwell's 
equations  in  a  homogeneous  source-free  region. 


TE  Fields  TM  Fields 


Ep  -  - 

(j?i3u/p) 

RP'  Z 

Zp 

m 

r'  p  z 

E*‘ 

juu 

R*  P  Z 

E* 

m 

(1/p)  R  P:  Z 

Ez  “ 

0 

Ez 

m 

B2  R  P  Z 

Hp  - 

r'  p  z' 

HP 

& 

(juc/p)  R  P*  Z 

H«- 

(1/p) 

R  P'  Z’ 

** 

m 

-jue  R1  PZ 

Hz  - 

82 

R  P  7. 

Hz 

- 

0 

The  tine  dependence  e^ut  is  understood.  R  is  a  function  of  p  only, 

P  is  a  function  of  4>  only,  and  Z  is  a  function  of  z  only.  Primes  indicate 
differentiation  with  respect  to  p,  ^  or  i,  The  functions  satisfy  the 
following  differential  equations: 

P  +  (S2  P2  -  m2)  R  -  0 

OP 

P"  -  -  m2  P 
Z"  -  -  h2  Z 


where 


,  2  •) 

8^  +  h  ■  id  vc,  and  8  and  h  are  constants. 


Some  solutions  of  these  differential  equations  are  listed  below. 
R(p)  *  Jo  (8p)  P($)  ■  cos  Z(z)  *  cos  hz 


No  <BP> 

sin 

sin  hz 

(Bp) 

ejo« 

ejhz 

Hi" 

(Bp) 

-jw> 

e 
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APPENDIX  II 


The  geodesic  differential  equations  for  an  arbitrary  surface  of 
revolution  are  developed  in  this  section  using  tenser  analysis.  Using 
the  geometry  illustrated  in  Fig.  56,  any  point  on  this  surface  is 
defined  by 

x(e,o)  =  R(e)  sin  e  cos  $ 

(106)  y(e,$j  =  R(e)  sin  e  sin  $ 
z(e,*)  =  R(e)  cos  e. 

The  unit  tangent  vectors  on  the  surface  are  given  by 
e-jfe,*)  =  [x  sin  e  +  r  cos  e^cos  $  +  y^jj|  sin 

(107)  2  cos  0  -  R  sin 
and 


A  A  A 

e2(<J>><p)  =  -x  sin  <j>  +  y  cos  <J>. 
The  metric  tensor  is  given  by 


Using  the  above  information,  the  Christoff el  symbols  are  given  by 

j2r 


r, 

11  if  * 


r1  = 
r22 


dR 

-  Rsin  e  (R  cos  e  +  ^  sin  e) 


(H) ♦  »2 


dR 

2  2  R  cos  0  +  cTe  S1*H  e 

r2i  =  r12  =  R  s>n -0 


rl  =  1  =  r2  2  =  Q 

r12  r21  rll  r22  U  • 


Substituting  these  results  into  the  geodesic  differential  equations^], 
one  finds  that  the  geodesics  paths  are  defined  by 


where  £  is  the  arclength  along  the  geodesic  path.  Note  that  these 
solutions  are  valid  only  for  a  surface  of  revolution  such  that  the 
z-axis  coincides  with  the  axis  of  revolution. 
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